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ABSTRACT 


A  recent  review  of  nucleer  ground ehock  test  data  revealed  oecillatory  notions 
of  undetermined  origin,  but  of  sufficient  magnitude  to  be  of  concern  In  the 
design  of  shock  isolation  systems  for  underground  protective  structures. 
Response  spectra  of  the  oscillations  exhibited  amplification  ratios  far 
exceeding  those  employed  In  most  groundshock  prediction  methods,  the  ratios 
being  functions  of  the  number  of  oscillations,  and  the  amplitude  and  period  of 
each  cycle.  As  the  source  of  the  oscillations  was  not  Identified,  there  was 
no  rational  basis  for  relating  these  waveform  parameters  to  such  basic  site 
and  weapon  conditions  as  yield,  range,  geologic  structure  end  properties  of 
the  medium.  In  this  Investigation,  two  possible  sources  of  the  oscillatory 
motion  are  examined.  First,  the  propagation  of  waves  In  a  stratified  site 
are  studied  and  their  directions  and  phase  reletlonships  estimated  by  use  cf 
wavefront  diagrams  and  time-distance  curves.  Second,  the  form  end  strength  of 
Rayleigh  waves  In  an  elastic,  homogeneous  half-space  w'uich  res  It  from  surface 
pressure  distributions  similar  to  those  generated  by  nucleer  bursts  were 
calculated.  In  both  cases,  oscillatory  phenomena  can  be  predicted  and  certain 
features  related  to  the  observed  oscillations.  However,  the  simple  approaches 
employed  in  this  analysis  will  not  yield  realistic  wave  strengths  end  thus, 
the  composite  vsveform  at  a  point  In  the  half-space  cannot  be  determined 
quentitatlvely. 
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i.  mrooDucTioN. 

A.  Background.  For  the  past  several  years ,  considerable  attention 
has  been  focused  on  various  aspects  of  the  design  of  shock  isolation 
systems  for  use  in  undergo ound  protective  structures.  Abetter  definition 
of  the  groundehock,  of  course,  has  been  the  goal  of  numerous  studies  of 
weapons  effects,  both  from  the  theoretical  approach  and  from  continued 
efforts  to  extract  more  detailed  Information  from  the  limited  number  of 
nuclear  test  data,  between  the  free  field  ground  motion  and  the  input  to 
the  shock  isolation  system  lies  a  capsule  or  other  enclosing  structure, 
and  its  dynamic  response  to  the  ground shock  defines  the  nature  and  strength 
of  the  shock  from  which  equipment  and  personnel  must  be  protected.  Here 
again,  both  theory  and  experiment  are  being  employed  to  evaluate  the  inter* 
actions  and  to  formulate  design  procedures. 

The  analysis  of  the  shock  isolation  system  Itself  as  an  elastically  sup* 
ported  rigid  body  having  six  degrees-of-freedon  has  been  the  subject  of 
Innumerable  studies;  although,  in  fact,  methods  for  treating  such  problems, 
both  linear  and  nonlinear  have  been  well  established  for  some  time.  In 
this  case  the  objective  has  bee  more  one  of  familiarising  the  designer 
with  available  techniques  than  of  developing  new  ones.  Similarly,  the 
feasibility  of  a  wide  variety  of  shock  isolation  devices,  generally  based 
on  adaptations  of  existing  principles,  has  been  Investigated  and  various 
systems  have  been  proposed  as  offering  unique  advantages. 

With  the  accumulation  of  a  siseable  fund  of  operating  and  teat  experience, 
and  here,  of  course,  "operating "  laplies  static  operation  and  not  exposure 

to  the  design  shock,  more  detailed  requirements  bars  been  encountered  and 

-  * 

reported  in  the  literature.  Such  considerations,  for  exwjtle,  include 
ccapensation  for  load  changes,  maintenance  techniques  which  can  be  employed 
without  sacrificing  hardness  level,  significance  of  msnufscturlng  tolerances 
on  performance,  evaluations  of  toe  magnitude  of  air  loads  acting  on  the 
structure  due  to  Its  motion  relative  u>  the  capsule,  effects  of  flexibility 
of  the  eupported  structure  and  the  transmission  of  high-frequency  accele¬ 
rations  through  the  amcemulcal  elements  of  tbs  isolator. 
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Concept  studies  of  advanced  facilities  have  also  contributed  to  the  recog¬ 
nition  of  problems  relating  to  the  design  of  current  as  veil  as  future 
•bock  Isolation  system  and  have  stimulated  the  investigation  of  exotic 
designs.  %  extrapolating  current  designs  to  very  large  or  very  hard 
installations!  weaknesses  in  the  designs  and  the  design  procedures  now 
•nplo^'?d  frequently  have  bee  one  apparent. 


Whan  each  of  these  problem  areas  Is  examined  in  detail  and  it  is  ascer¬ 
tained  whether  engineering  information  or  engineering. knowledge  is  needed , 

It  appears  to  be  clear  that  almost  all  of  the  problems  fall  in  the  former 
category.  It  is  not  lulled  that  the  design  of  shock  isolation  systems 
for  this  application  has  been  reduced  to  pure  routine  and  that  no  further 
developmental  or  research  effort  Is  warranted.  It  is  believed ,  however* 
that  if  performance  criteria  are  established  completely  and  quantitatively 
prior  to  design*  and  if  a  development  program  la  conducted  in  accordance 
with  the  high  standards  practiced  in  equally  important  engineering  under¬ 
takings*  most  of  the  problems  noted  above  could  be  resolved  with  a  high 
level  of  confidence. 

The  outstanding  exception  is  with  those  elements  of  tae  isolation  system 
which  are  sensitive  to  aspects  of  the  input  which  cannot  be  evaluated  with 
the  needed  accuracy  by  any  known  weapons  effects  prediction  method.  For 
example*  certain  features  of  the  tins  history  of  shock  must  be  known 
in  order  to  determine  the  dynamic  responses  of  nonlinear  or  r he oil near 
elements  or  systems.  Sven  with  coupled  or  hig;Jy  damped  linear  systems* 
accurate  calculation  of  tbs  peak  responses  requires  some  information  on  the 
nature  of  the  waveform*  and*  to  evaluate  the  frequency  content  of  the  output* 
the  Input  disturbance  must  be  known  as  a  function  of  time.  While  a  number 
of  structural  elements  and  a  few  shock  isolation  systems  may  be  ereroxi- 
mated  with  acceptable  accuracy  as  mingle  degree-of -freedom,  lightly  duped 
oscillators*  a  mast  significant  percentage  csanc*  lima*  the  identification 
of  input  waveform  parameters  which  oaa  Influence  the  response*  of  systems 
typical  of  those  employed  in  underground  protective  *ructures  appears  to 
be  the  moat  important  single  research  problem  remaining  to  be  resolved. 
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It  is  possible  through  the  judicious  selection  of  suspension  kinematics 
and  isolator  characteristics  to  minimise  this  dependence  of  isolation 
system  response  on  the  less  predictable  features  of  the  input  waveform. 
Indeed ,  to  ensure  that  the  highest  level  of  confidence  in  the  ability  of 
the  system  to  perform  as  expected  is  achieved/  practical  designs  should 
begin  vith  insensitivity  as  one  of  the  stated  design  objectives.  By 
assigning  levels  of  confidence  to  each  of  the  predicted  parameters  of  the 
shock  input  and  by  relating  the  various  responses  or  the  system  to  these 
parameters/  the  direction  in  which  optimization  should  proceed  becomes 
clearly  evident. 

However/  confidence  level  is  not  the  only  criterion  of  performance,  and 
other  design  goals  render  it  rarely  possible  to  achieve  the  decree  of 
uncoupling  of  response  from  input  which  the  uncertainty  surrounding  the 
details  of  the  input  would  appear  to  warrant.  Thus  while  a  strong  effort 
should  be  aade  in  every  case  to  reduce  the  sensitivity  of  the  system  to  a 
minimum  consistent  with  the  other  performance  objectives/  the  design  of 
most  practical  systems  will  require  more  information  concerning  the  input 
than  tnat  afforded  by  current  groundshock  prediction  techniques 

Until  recently,  the  only  information 'describing  the  nature  of  the  shock 
environment  which  was  specified  as  a  basis  for  the  Assign  of  shook  iso* 

1st Ira  systems  was  a  frequency  spectrum  of  the  envelope  of  the  response* 
of  linear/  undamped  systems.  Sheet  spectra  were  formulated  by  groundshock 
prediction  methods  based  on  tests  conducted  at  the  Arvada  fast  Sits  and  the 
Sol  vet ok  Proving  Ground.  Both  waveforms  and  spectra  mure  obtained  during 
the  teats.  However  only  gross  characteristic  features  of  tbs  waveforms 
coul*  be  | identified;  thus,  except  far  a  Halted  number  of  paymasters,  the 
prediction  methods  were  directed  toward  the  ayntbesls  of  response  spectra 
rather  then  of  time  histories,  further,  the  geologic  formations  at  the 
teat  sites  ears  not  typical  of  tho-e  at  which  moat  underground  protective 
structures  are  constructed,  nor  ware  the  weapon  sizes  of  the  igrtltwde  of 
therm  now  being  conslAsrad.  In  addition/  insufficient  site  end  weapon 
parameters  could  he  vari  A  during  the  testa  to  evaluate  independently  the 


•ffieti  of  eaeb.  Tbarefcs re,  to  extrapolate  these  data  to  otter  site  and 

conditions  of  ii  merest,  sljqple  linear >  one-d.imensional  relationships 

used  extensively. 


Ai  a  Beans  for  estimating  the  gross  strength  of  the  shock  which  might  be 
SBgietaft  to  ftfleur  under  pacified  conditions ,  these  prediction  methods  have 
felia  of  tltal  importance  to  the  protective  structure  program.  Indeed  with- 
inf  thilij  it  is  unlikely  that  the  design  of  structures  of  significant 
bsrdrmrs  would  tana  been  attempted,  thus  these  prediction  techniques , 
orala  as  they  any  be*  have  played  an  invaluable  role  in  the  construction 
of  all  bardsaad  facilities  in  the  United  States. 


Boaavar*  to  utilise  the  limited  information  available  on  the  nature  of  the 
graundahoek*  early  analyses  of  shock  isolation  systems  and  other  structural 
components  of  protective  facilities  esployed,  of  necessity,  many  siapUfying 
assumptions.  In  most  cases  the  systems  were  regarded  as  being  li  near  and 
vbtro  histories  of  the  input  were  essential,  any  pulse  whose  response 
•paetrw  matched  or  exceeded  the  design  spectrum  in  the  frequency  range  of 
Interest  was  accepted  aa  a  suitable  representation.  However  as  the  intended 
level  of  each  new  facility  increased,  the  requirements  for  an  ever 
greater  reliability  and  a  more  nearly  optimum  design  rose  correspondingly, 
dpelyaea  more  rigorous  and  now  have  been  extended  to  many  elements 

formerly  treated  only  by  groea  approximation.  In  almost  every  case  the 
greeter  emphasis  on  rigor  In  design  analysis  has  demanded  a  more  accurate 
«ai  detailed  definitive  of  the  shock  environment. 

fhr  the  moat-  part,  currant  weapons  affects  prediction  methods  consider  two 
basic  grotaaRahock  pulse  shapes ,  the  airb last-induced  pulse  and  the  direct- 
in'uced  pulse*  Procedures  are  givan  for  calculating  the  posh  motions  for 
amah  type  of  pulse  and,  together  with  suitable  amplification  ratios,  these 
yaak  values  in  used  as  a  basis  for  constructing  response  spectra.  Some 
j|i  “tmrna  is  provided  for  eetlmatlx*  certain  other  parameters  of  their  wave¬ 
form  tag  am  acceleration  rise  time,  positive  velocity  duration,  and  total 
duration,  the  confidence  level  in  the  predictions  of  these  latter 
p— bars  bomavar  is  wwaliy  considered  to  be  far  leas  than  those  of  the 
yaak  motions. 
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To  account  for  the  effects  of  all  other  ground  notions  generated  by  the 
nuclear  burst,  and  in  particular  of  the  reflections  and  refractions  from 
the  various  layers  of  stratified  media,  the  usual  practice  is  to  increase 
the  peak  airblast-induced  motions  by  a  factor  vhich  is  related  to  the  range 
of  the  point  of  Interest  and  the  range  at  which  outrunning  first  occurs. 
While  LI 1-  method  may  yield  acceptable  values  for  the  peak  ground  motions, 
it  given  no  indication  of  the  detailed  nature  of  the  waveform  other  than 
to  describe  it  as  a  "random  oscillation."  As  a  great  many  facilities  of 
interest  are  constructed  at  highly  stratified  sites  the  problem  of  defining 
time  histories  of  the  ground shock  is  an  important  one. 

In  the  design  of  the  shock  isolation  system  fen'  the  launch  Control  Center 
of  the  KLmitemsn  Weapon  System,  the  ccoAined  requirements  for  attenuation 
and  rattlespace  were  impossible  to  satisfy  with  a  linear  system,  aixJ  as 
only  a  response  spectrum  was  provided  to  define  the  shock,  the  engineering 
contractor  was  faced  with  the  responsibility  for  synthesizing  waveforms 
which  could  be  esployea  to  verify  the  performance  of  the  nonlinear  system 
he  had  selected.  Considerable  discussion  then  centered  about  the  validity 
of  the  waveforms  be  had  constructed  and  there  nlway  remained  some  question 
that  the  most  critical,  conditions  had  been  represented. 

Partially  cm  the  basis  of  this  experience,  the  Air  Force  Special  Weapons 
Center  in  19Sl  initiated  several  studies  directed  toward  the  formulation  of 
design  methods  applicable  to  nonlinear  shock  Isolation  systems  for  use  in 
underground  protective  structures.  In  one  of  these  (l),  all  available 
waveform  data  from  nuclear  tests  were  re -examined  with  the  objective  of 
Identifying  all  characteristic  features  of  the  waves  which  might  be  of 
significance  to  the  responses  of  low  frequency  systems.  It  was  found  that 
when  the  contribution  of  the  direct  airblast-induced  motion  was  removed 
from  certain  test  records,  there  remained  a  low  frequency  oscillation  of 
sufficient  strength  to  be  of  vital  concern  In  the  design  of  linear  as  well 
as  nonlinear  shock  isolation  systems.  Although  the  oscillation  could  not 
be  traced  directly  to  its  source,  it  was  suggested  that  it  resulted  from 
refractions  returning  to  the  surface  from  underlying  layer*  having  higher 
seismic  velocities.  This  hypothesis  appeared  to  correlate  well  with  the 
fact  that  In  the  test  data,  the  oscillation  wax  most  apparent  Is  the 
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txans-seismic  range  where  the  refract  lone  vould  be  expected  to  contribute  a 
large  percentage  to  the  total  ground  motion. 

To  explore  further  the  possible  source  of  such  an  oscillation  the  present 
investigation  was  undertaken.  Other  research  projects  now  in  progress  are 
Studying  the  general  problem  of  the  transmission  of  waves  through  homoge¬ 
neous  anf  layered,  elastic  and  inelastic  media.  This  effort,  however,  was 
concerned  solely  with  those  aspects  of  groundshoclc  phenomena  which  might 
shed  some  light  on  the  source  and  strength  of  the  observed  oscillations. 

The  problem  of  identifying  possible  sources  of  the  waves  was  approached  by 
two  distinctly  different  methods.  First,  the  wavefront  diagram  was  employed 
to  determine  whether  or  not  the  oscillations  appeared  to  be  simple  combi¬ 
nations  of  airb last -induced ,  direct -induced,  reflected  ^nd  refracted  waves 
Inpinging  at  a  point  in  some  sequence  defined  by  the  geology  of  the  site. 

Hie  wavefront  diagram  is  a  graphical  representation  in  a  vertical  plane  of 
a  wavefront  as  it  propagates  outward  from  a  point  source,  and  although 
Inelastic  media  and  arbitrary  stratification  of  the  site  can  be  considered, 
it  is  evident  that  the  method  is  strictly  limited  in  its  ability  to  portray 
faithfully  all  wave  phenomena  emanating  from  a  nuclear  explosion.  Nonethe¬ 
less,  as  noted  earlier,  the  objective  here  was  only  to  establish  the  phase 

relationships  and  if  possible  the  relative  strengths  of  the  four  types  of 

« 

waves /and,  in  this  respect,  the  wavefront  diagram  appeared  to  offer  some 
promise. 

In  the* second  approach,  the  contribution  to  the  oscillatory  motion  of  surface 
waves  were  investigated  theoretically.  Again,  however,  only  forms  and  the 
relative  magnitudes  of  such  waves  need  be  ascertained  to  class  them  as  likely 
relatives  of  those  observed  in  the  test  data.  Thun,  this  effort  also  was 
conducted  without  the  degree  of  detail  considered  necessary  in  a  more 
rigorous  analysis  of  groundshock  wave  propagation.  In  particular,  the  form 
and  strengths  of  Rayleigh  waves  generated  in  an  elastic,  homogeneous  half- 
epace  by  a  surface  pressure  distribution  similar  to  that  resulting  from  a 
nuclear  hurst  were  calculated  and  their  phasing  with  other  waves  emanating 
from  the  source  estimated.  They  were  then  cospared  with  the  observed 
oscillatory  motions  and  points  of  similarity  noted. 
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B.  Observed  Oscillatory  Groundshock.  Sauer  (l*  25)  was  the  first  to 
attempt  to  synthesize  the  general  characteristics  of  the  strong  velocity 
oscillations  observed  in  several  groundshock  waveforms  integrated  from 
acceleration  records  obtained  during  nuclear  weapon  tests.  While  high- 
frequency  oscillatory  phenomena*  of  course*  are  frequently  found  in 
acceleration  traces  of  this  type,  they  cam  usually  be  expected  to  have 
little  influence  on  the  ..velocity  curve.  In  this  case*  the  relatively  low 
frequency  of  the  velocity  oscillations  together  with  their  significant 
strength  and  apparent  regularity  of  form  lead  to  the  hypothesis  that  they 
resulted  from  some  undetermined  but  not  unusual  interaction  of  the  various 
groundshock  components. 

The  oscillations  were  observed  in  data  taken  both  at  the  Enivetok  Proving 
Ground  (EFG)  and  at  the  Nevada  Test  Site  (NTS).  To  isolate  the  oscillation, 
the  essentially  regular  one-sided  velocity  waveform  commonly  associated 
with  the  superseismic  airblast-induced  groundshock*  designated  as  a  Type  I 
wave  and  shown  in  Figure  1*  was  subtracted  from  outrunning  Operation  Tumbler 
data  and  compared  with  data  from  Xoa  12  where  the  Type  I  was*  presumably, 
filtered  out  by  the  ground.  The  relative  amplitudes  of  the  various  peaks 
were  averaged  over  these  data  as  were  the  relative  tine  durations  of  the 
cycles.  More  weight  was  given  to  the  Koa  data  since  the  yield  for  this  shot 
was  in  the  megaton  range.  After  the  first  tentative  velocity  curve  had  been 
constructed  it  was  integrated  to  give  the  displacement  and  adjustments  were 
made  in  the  later  portions  of  the  velocity  curve  so  that  the  residual  dis¬ 
placement  would  be  zero. 

The  resulting  oecillation,  designated  as  a  Type  II  wave*  is  shewn  in 
Figure  2  and  is  compared  in  Figures  5  and  4  with  several  traces  obtained  at 
Koa.  In  some  cases  the  normalized  curve  fits  the  data  very  well;  in  others* 
the  fit  is  not  so  good,  but  the  normalized  curve  still  has  the  general  shape 
of  the  measured  motions.  Sauer  suggests  that  the  oscillation  is  related  in 
some  fashion  to  the  outrunning  coalition  and  illustrates  the  development  of 
outrunning  waveforms  obtained  during  Operation  Tumbler*  Shot  No.  1*  by 
appropriate  combinations  of  T/P®  I  and  typ®  II  waveform*  (25). 

In  reviewing  the  test  data,  Sauer  noted  that  no  characteristic  features 
could  he  identified  in  the  acceleration  records.  It  seems  reasonable  hew- 
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ever  that  if  the  Ifype  II  oscillation  is  the  interaction  of  several  waves 
reaching  the  point  of  interest  by  different  routes  and  at  different  times , 
the  arrival  of  each  wave  should  be  indicated  by  a  fairly  abrupt  change  in 
acceleration .  It  is  possible  t ' though >  that  the  number  of  waves  resulting 
frm  the |  disturbance  in  a  highly  stratified  >  nonlinear  medium  was  so  large 
t&at  their  identity  was  lost. 
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Figure  1.  Type  I  (Super seismic  Airblast)  Normalized  Vertical 
Particle  Velocity  and  Displacement  Waveforms. 


Figure  2.  Type  II  (Outrunning  Ground  Motion)  Horstlized  Vertical 
Particle  Velocity  and  Displacement  Wavs  *oru. 
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Figure  3.  Comparison  of  the  type  II  Vertical  Velocity 
Waveform  with  Data  from  Shot  Cactus,  Gage 
2V30,  630  Foot  Range,  30  Foot  Depth. 


Figure  k.  Co^arison  of  the  type  n  Vertical  Velocity 
Waveform  with  Data  from  Shot  Koa,  Ga<p 
12V1C0A,  3#lW»  Foot  Range,  100  Foot  lepth. 

Nonetheless  it  was  believed  to  be  sufficiently  promising  1;o  investigate  the 
possibility  of  constructing  a  reasonably  close  approximation  of  the  shape 
of  the  type  n  component  by  employing  information  obtained  from  the  simple 
wavefront  and  time -distance  diagrams  together  with  the  few  data  provldad  by 
the  groundshock  prediction  methods.  Although  many  amplifications  would  ha 
required.  It  was  hoped  that  the  most  Jjpcrtant  phenomena  could  be  considered 
and  that  the  resulting  waveform  would  be  a  significant  isprovement  over  the 
arbitrarily  selected  shapes  previously  alloyed  in  shock  isolation  system 
analyses. 
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A*  Background.  Calculation  of  the  transmission  of  energy  from  an 
air  or  surface  burst  of  a  nuclear  weapon  to  a  given  point  on  or  In  the 
ground  la  an  exceedingly  difficult  task*  even  for  sites  whose  geologic 
structure  is  relatively  uniform.  While  several  investigators  have  form¬ 
ulated  coaguter  programs  to  calculate  the  ground  motions  In  inelastic, 
layered  media,  to  date  no  computer  solutions  have  provided  time  histories 
of  the  ground  motion  suitable  for  design  purposes. 

On  the  other  hand,  accepted  weapons  effects  prediction  methods  yield 
little  more  than  naxlmum  motions  of  the  ground  resulting  from  a  single 
wave.  Slaqple  pulse  type  waveforms  are  suggested  for  airblaat- induced 
and  direct -induced  waves,  hut  no  guidance  is  given  for  predicting  wave¬ 
forms  which  might  result  hy  the  impingement  at  a  point  of  several  waves 
of  different  magnitudes,  from  different  directions  and  at  phase  intervals 
related  to  the  specific  site  under  consideration.  Almost  all  of  these 
works  have  been  directed  toward  the  synthesis  of  response  spectra  of  the 
groundabock,  not  a  waveform,  and  thus  their  direct  applicability  to  shock 
Isolation  system  design  has  been  limited  to  linear  systems  whose  equations 
of  motion  contain  only  constant  coefficients. 

As  a  means  for  better  visualising  the  modes  of  propagation  of  the  shock 
through  the  ground,  the  df signer  may  esgloy  a  wavefront  diagram.  ‘The 
wavefront  diagram  is  singly  a  representation  of  a  vertical  plans  passing 
through  the  point  of  weapon  burst  and  the  site  of  interest  on  which  the 
position  of  the  wavefront  Is  shewn  at  discrete  time  intervals,  Any  number 
of  layers  can  be  shown  and  oriented  in  whatever  fashion  the  designer 
believes  to  be  most  representative  of  the  actual  site,  the  degree  of 
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detail  that  can  be  included  in  the  wavefront  diagram  is  limited  prirau  ily 
by  the  patience  and  resourcefulness  of  the  designer,  but  in  most  instances 
the  assumptions  of  elastic,  homogeneous ,  isotropic  layers  with  sharp 
interfaces  have  been  employed. 

Data  obtainable  from  a  wavefront  diagram  are  primarily  the  arrival  time  of 
a  wave  at  the  point  of  interest,  the  path  by  which  the  wave  reached  the 
point  and  the  initial  direction  of  propegs, t ion  of  the  wavefront.  If  the 
assumption  of  a  linear  medium  either  with  or  without  layering  is  reasonably 
valid,  these  data  can  be  determined  for  all  waves  generated  directly  by 
the  source  or  by  the  me  ring  alrblast  wavefront.  In  most  cases  of  interest 
however  the  medium  is  not  linear  even  at  ranges  extending  beyond  the 
rupture  zone  and  thus  the  velocity  of  propagation  of  each  wave  at  each 
instant  is  dependent  on  the  local  time-stress  history  of  the  medium.  Sven 
if  the  nonlinear  stress-strain  characteristic  of  the  medium  is  known,  the 
strength  of  the  waves  cannot  be  determined  from  the  wavefront  diagram  and 
thus  the  propagation  velocity  of  waves  of  finite  strength  can  only  be 
approximated. 

3he  first  objective  of  this  phase  of  the  work  then  was  to  coe^are  computed 
arrival  times  of  various  waves  for  site  and  weapon  conditions  corresponding 
+o  those  prevailing  at  selected  nuclear  tests  with  the  measured  arrival 
times  in  order  to  determine  the  degree  of  error  involved  in  the  assumption 
0  of  a  linear  medium.  For  a  solution  of  the  second  problem,  that  is  the 

estimation  of  the  time  of  arrival  and  the  magnitude  of  the  peak  motions  of 
each  wave,  the  wavefront  and  the  ti^-distance  diagrams,  of  course,  provide 
little  useful  data. 

B.  Tbs  Use  of  the  Ihvefront  Diagram  in  Ground  Motion  Problems. 

The  basis  and  procedure  for  constructing  wavefront  diagrams  and  accom¬ 
panying  tins -distance  curves  of  tbs  refracted  and  reflected  waves  and 
nethods  of  using  these  diagrams  with  the  alrblast  curve  to  compute  the 

11 
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time  of  arrival  of  the  firrt  wave  are  presented,  in  this  section.  Ifce 
Methods  are  then  applied  to  test  shots  made  during  Operation  Tumbler  and 
Shot  Priscilla  of  Operation  Plusfcbob.  Seismic  velocities  measured  during  a 
conventional  seismic  survey  are  compared  with  the  velocities  measured  during 
the  test  shots  for  the  purpose  of  determining  the  effect  of  nonlinearity  of 
the  ooll  and  establishing  a  basis  for  computing  the  time  of  arrival  of  waves 
subsequent  to  the  first.  Finally,  the  distribution  of  energy  at  the  inter¬ 
face  of  two  elastic  layers  is  examined  to  determine  the  relative  proportions 
of  energy  carried  by  the  reflected  and  refracted  waves. 

(l)  The  Wavefront  Diagrea.  Wavefront  diagrams  were  developed 
as  a  tool  in  seismic  prospecting  to  obtain  a  better  "feel"  tar  the 
interpretation  of  the  time -distance  graphs  obtained  from  refraction 
seismic  surveys.  Examples  of  the  use  of  wavefront  diagrams  appeared  in  the 
literature  as  early  as  1950,  (2).  Wavefront  diagrams  provide  a  visual 
method  for  evaluating  both  qualitatively  and,  with  regard  to  time  and 
distance,  quantitatively,  the  nature  of  the  motion  of  the  ground  when  a 
source  of  energy  \s  applied  at  a  point  on  the  surface. 

(A)  Baals.  A  typical  wavefront  diagram  ie  presented  in 
Figure  %  The  diagram  is  singly  a  composite,  planar,  vertical  cross* 
section  of  the  fraud  showing  the  positions  of  the  wavefront  emanating 
from  a  point  source  at  specific  tlae  intervals.  At  any  given  tine  the 
wavefront  is  a  curved  surface  passing  through  the  most  advanced  positions 
reached  by  the  disturbance  at  that  time. 

lb  construct  a  wavefront  diagram,  it  is  convenient,  although  not  essential, 
to  assume  that  the  wave  velocity  in  each  layer  or  term  is  constant.  To 
coagwte  arrival  tines  no  assumption  is  necessary  regarding  the  character 
of  the  disturbance.  The  boundaries  of  the  various  layers  or  tone*  nay  be 
irregular,  or  oriented  in  any  direction,  and  the  wave  propagation  velocity 
can  be  allowed  to  change  in  small  or  large  steps,  and  in  either  or  both 
the  vertical  or  horizontal  directions. 
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For  any  arrangement  of  layers  with  different  wave  velocities,  the  entire 
net  can  be  constructed  using  Huygen's  principle:  every  point  on  a  wave- 
front  acts  as  though  it  were  itself  a  center  of  disturbance,  sending  out 
wavelets  of  its  own,  always  away  frees  the  source,  the  collective  effect 
of  which  constitutes  a  new  wavefront.  For  coeplex  profiles,  it  Is  usually 
easier  to  construct  the  wavefront  diagram  using  Just  the  basic  principle  of 
Huy gen;  however,  for  regular  profiles,  construction  can  be  speeded  by 
certain  other  aids.  A  detailed  step-by-step  example  of  the  construction 
of  a  wavefront  diagram  is  given  in  Appendix  A. 


(B)  Character.  All  of  the  wavefronts  shown  in  Figure  5 
are  not  just  simple  spherical  waves  emanating  from  a  point,  but  amj  be 
composed  of  several  wavefronts  of  different  origin  and  history.  For 
example,  the  wavefront  at  a  time  equaling  3.0  seconds  consists  of  a 
combination  of  four  simple  wavefronts.  The  section  of  the  front  traveling 
in  Layer  1  is  the  "first  underlayer  nave"  which  has  arrived  at  that 
position  by  being  refracted  along  the  top  of  layer  2  (first  underlayer) 
and  back  into  layer  1.  The  part  of  the  front  between  the  top  of  layer  2 
and  the  curve  BGS  is  the  refracted  wave  traveling  through  the  body  of 
layer  2.  Bie  remaining  part  of  the  front  in  layer  2  has  been  refracted 
into  Layer  3  and  back  into  layer  2.  This  is  called  the  "second  underlayer 
wave. *  The  part  of  the  front  in  Layer  3  i*  the  refracted  wave  traveling 
through  the  body  of  that  layer. 

Due  to  the  higher  velocity  of  layer  2  compared  to  Layer  1,  the  area  to  the 
right  of  Point  A  in  Figure  5  is  first  reached  and  set  in  motion  by  the 
refracted  wave  traveling  in  and  along  the  top  of  layer  2,  not  by  the  direct 
wave  from  Foist  $.  The  eagle  which  the  underlayer  wave  nahos  with  the 
boundary  at  which  it  originate#  is  independent  of  the  position  of  the 
original  s antes  or  the  formation  thickness  end  depends  only  on  the  velocity 
ratio  of  the  two  formation*  adjoining  the  boned  ary. 
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If  layer  2  had  a  velocity  lower  than  that  of  layer  1*  the  direct  waves  in 
Layer  1  would  be  refracted  downward  into  layer  2  and  no  "underlayer  wave" 
would  be  generated.  Only  for  those  boundaries  at  which  the  velocity 
increases  vith  depth  will  underlayer  waves  be  generated. 

Along  the  top  of  Layer  2,  the  wavefronts  to  the  right  of  Point  A  appear  to 
have  originated  at  Point  A.  The  ease  condition  exists  to  the  right  of 
Point  B  in  layer  J.  At  both  Points  A  and  B,  total  reflection  from  a 
source  at  S  would  occur,  according  to  Snell's  Lev. 

(C)  Shell's  law.  The  refraction  of  waves  at  the  layer 
interfaces  is  governed  by  Shell's  law.  This  law  was  discovered  export- 
aentally,  but  it  actually  is  consequence  of  Peraet's  principle  which 
states  that  rays  travel  along  such  lines  that  the  optical  distance 
between  any  two  points  of  the  rays  Is  a  Minim.  If  the  given  point  is 
in  a  MsdiuM  characterised  by  a  velocity  different  fron  the  Msdlun 
containing  the  second  point,  then  the  rey  path  between  these  two  pedate 
will  not  be  a  straight  line. 


Figure  6.  Angles  of  Incidence  aad  Refraction 


I 


Snell's  Ira  of  refraction  stater  that  refraction  or  bending  of  the  ray 
fgfifc  oeeors  when  cresting  a  bouw’ary  between  I aye**  of  different  velocities 
erafc  that*  raftering  to  figure  6* 


Ubeat  e  rgr  in  e  lower  velocity  radius  strikes  ea  interface  vlth  a  higher 
velocity  radius*  there  ia  s  certain  critical  angle  of  incidence  deslg- 
-  rated  as  1.*  for  which  the  eagle  of  refraction  r  is  90  degrees  and  the 
refracted  save  le  parallel  to  Idee  surface  of  discontinuity.  Since,  in 
this  ease*  sia  r  «  1*  it  Is  evident  free  Snell's  lees  that 


sin  i 
c 


(3) 


fbr  any  eagle  of  Incidence  greeter  thee  le*  there  can  be  no  refracted  ray 
la  the  sac  and  endive and*  therefore,  po  penetration  Into  that  nwliun,  all 
the  isergy  incident  r-t  steb  »  angle  being  reflected.  Thus,  for  incident 
eagles  greater  *haa  the  critical  angle  there  is  total  reflection. 


Point*  A  and  B*  In  Mger*  5,  represent  the  location  of  a  ray  path  frost 
MUt  S  at  the  critleel  eagle*  The  save  refractions  which  occur  to  the 
left  of  ^eftaffca  A  sad  1  do  aot  return  to  the  surface.  The  fact  that  Points 
A  sad  »  era  bs  loe&tad  oqr  Saell's  las  and  represent  the  aialnua  tioe  paths 
is  of  eooelderable  assistance  ia  tht  graphics!  construction  of  e  ravefront 
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(D)  Coiacldent-Tine  Curves.  Curves  represented  hi  points 
AK  and  BSD  of  Figure  5  are  for  and  wherever  two  different  wavefront 
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Aysteas  Intersect.  Together  with  the  layer  bouodAri.es ,  these  curves  divide 
the  section  Into  tones#  each  tone  confining  wavefronts  of  *  c canoe  origin. 
Bsr  construction#  ell  points  an  such  &  curve  are  equally  distant  in  tine 
iron  the  source  along  tvo  different  paths.  Therefore#  such  lines  are  called 
coincident -tine  curves. 


(K)  later  Arrivals.  The  wavefronts  shewn  in  Figure  5  are 
only  those  arriving  first  at  each  point.  The  direct  wave  shown  as  the 
first  arrival  to  the  left  of  line  AFC  in  Figure  5  will  continue  arriving 
later  at  points  C#  X,  L#  D,  Z#  etc.  However,  this  wave  will  reach  these 
points  after  the  arrival  of  eaves  traveling  hy  other  pat)  The  saae  is 
true  with  the  first  underlayer  save  where  it  Is  no  longer  a  first  arrival 
to  the  light  line  of  HD.  Deflections#  which  are  never  a  first  arrival# 
will  also  as  rive  later  at  all  paints  along  the  ground  surface.  In  addi¬ 
tion#  the  reflection*  will  occur  at  the  interface  between  layer  1  and 
layer  2,  and  also  between  layer  2  and  layer 

Since  Shear  wares  travel  at  a  velocity  of  about  one-half  that  of 
longitudinal  waves#  they  do  not  appear  in  a  first  arrival  wavefront 
dlagraa.  Shear  waves  will  be  generated  ever  a  longitudinal  wave 
strikes  a  discontinuity  where  there  la  a  change  of  elastic  properties. 

In  fact#  four  mm  waves  will  he  set  up  in  this  situation.  Far  exaaplc# 
Figure  7  shows  an  Incident  longitudinal  wave  F  and  the  resulting  reflected 
longitudinal  wave  ff,  >  reflected  transverse  wave  IS.  #  refracted  loagita- 
<U«X  «t. 

#  to  the  rays  (—*——)  Indicate  longitudinal,  waves}  those  perpendicular 

to  the  nqr*  (  ■■  )  indicate  transferee  waves. )  A  &iaHnr  sat  of  four 

waves  (hut  at  different  angles)  would  be  set  35  if  the  lncidsnrt  wave  ware 
transverse. 
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Figure  7.  Reflected  and  Refracted  Transverse  and  Longi¬ 
tudinal  Wives  due  to  an  Incident  Longitudinal 
Wave. 

The  division  of  energy  between  the  different  waves  depends  on  the  velo¬ 
cities  and  the  densities  of  the  two  media  (i.®.,  Vlpl  and  V^g)  and  oh  the 
angle  between  the  incident  wave  and  the  surface  of  discontinuity.  The 
energy  division  between  the  various  waves  will  be  discussed  in  a  later 
section.  In  general,  for  simple  reflection,  the  proportion  of  reflected 
energy  will  be  greater  for  greater  differences  in  the  elastic  properties 
on  the  two  sides  of  the  discontinuity. 


In  seismic  prospecting,  concern  is  always  with  the  fastest  waves  and, 
therefore,  with  incident  longitudinal  waves  P,  reflected  longitudinal 
waves  PFj^,  sud  refracted  longitudinal  waves  PPg.  At  the  present  st  *ge  of 
develop  sent  of  seismic  prospecting,  little  use  is  made  of  the  shear  waves. 
Occasionally,  however,  special  purpose  surveys  are  run  to  measure  the  shear 
wave  velocity  to  provide  data  to  calculate  the  elastic  properties  of  a 
material,  such  as  Poisson's  ratio,  (3). 
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The  later  arrivals  can  be  shown  on  a  wavefront  diagram,  but  their 
presentation  beconos  very  cumbersome.  They  are  sore  easily  presented  on 
time-distance  curves,  which  will  be  discussed  in  the  following  section. 

An  example  of  a  reflection  wavefront  diagram  is  ohcnm  in  Figure  58. 

(2)  Tiae -Distance  Curves.  Time -distance  curves  nay  be 
constructed  frou  data  obtained  frees  a  seismic  survey.  In  general,  a 
particular  survey  will  consist  of  measurements  of  either  the  refracted  and 
direct  waves  or  the  reflected  waves.  The  time  of  arrival  of  the  wave  is 
determine  from  inspection  end  correlation  of  the  seismogram  from  each  • 
geophone,  or  record  of  motion  vs.  tins.  The  location  of  3&sh  gffophooe  in 
relation  to  the  disturbance  or  "shot"  is  known.  Th»  most  common  type  of 
seismic  survey  used  for  site  evaluation  for  civil  engineering  purposes  is 
the  refraction  survey;  however,  the  reflection  survey  is  the  most  cannon  in 
petroleum  or  other  deep  studies.  It  will  be  shown  later  that  a  reflection 
survey  would  also  be  of  significant  value  in  the  evaluation  of  ground 
motions  due  to  nuclear  -- 

The  value  in  this  application  of  the  Ums-dista&ce  curves  lies  in  the 
ease  with  which  the  sequence  of  arrival  times,  including  the  first  and 
later  arrivals,  at  any  point  of  Interest  can  be  identified. 

The  direction  of  motion  at  the  point  of  interest  cannot  be  determined 
from  the  tire -distance  curves  and  requites  either  the  construction  of  a 
wavefront  diagram  or  a  construction  for  the  particular  wave  in  question. 

For  Interest,  the  tiaa -distance  curve  derived  from  Figure  5  is  presented 
in  Figure  8.  Additional  time-distance  curves  for  the  test  sites  at 
Frenchman  Flats  uni  T-7  (fucca  Fiats)  in  the  Bevada  Proving  Grounds  are 
Shown  in  Figure  2b  and  Figures  26  through  28. 
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(A)  Buis.  Figure  9  represents  a  subsurface  section  con¬ 
sisting  of  two  horisootal  layers  which  are  assumed  to  be  isotropic  and 
newly  homogeneous,  hawing  nominal  longitudinal  wave  velocities  V1  and  v£ 
in  the  upper  and  lower  lasers  respectively.  Assume  also  that  velocities 
Vg  increase  subtly  with  depth,  but  are  sufficiently  different  so 
that  thaze  is  an  abrupt  discontinuity  in  seismic  velocity  at  the  interface 
Point  0  represents  the  source  of  disturbance  and  S>  a  seismometer 
station  located  at  distance  x  from  0. 

Following  the  explosion  at  0,  surface  waves  and  three  -types  of  longitudi¬ 
nal  wwmrs,  direct?  reflected,  and  refracted,  are  received  at  3,  provided 
that  x  is  sufficiently  great.  The  direct  wave  travels  along  the  surface, 
following  path  OS;  the  reflected  wave  follows  the  path  CBS;  and  the 
refracted  wave  follows  the  path  QACS.  The  time -distance  curves  for  the 
three  longitudinal  waves  may  be  determined  as  follows. 


The  length  of  the  curved  path  06  is  approximately  equal  to  x,  the  slight 
curvature  being  due  solely  to  the  increase  of  velocity  with  depth.  Hence, 
the  time,  t,  required  by  the  direct  wave  to  traverse  the  distance  OS  **  x 


t  »  and  x  -  .tV^  (4) 

The  position  of  the  wavefront  on  the  time-distance  curve  (Figure  1.0 ) 
therefore  is  a  itralght  line  which  passes  through  the  origin  and  has  a 
slopn  of  magnitude  l/V^ 

From  the  geometry  of  Figure  9,  it  is  seen  that  the  distance  from  the  point 
point  0  to  the  point  of  reflection  B  is 


00  -  v®2  +  *2  -  I V^"1 
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tlhen  the  Interface  titwia  layer  X  and  Layer  2  is  parallel  to  the  surface, 
a  coalition  of  eymtxy  exists  and  the  distance  S®  »  B8.  The  total 
distance  traveled  by  the  reflected  save  is  then  OB  +  BS. 

s*b.  V** ♦  «>8 

Since  the  distance  (®  ♦  S)  is  also  equal  to  tV1>  the  tine  for  the 
reflected  wave  to  travel  distance  x  is. 


This  is  the  equation  for  a  rectangular  hyperbola,  as  plotted  in  Figure  1C. 

The  tine -distance  curve  of  the  vave  Which  is  refracted  near  t*e  critical 
angle,  and  therefore  travels  in  the  lover  layer  along  a  peth  which  is 
approximately  parallel  to  the  boundary,  is  obtained  by  adding  the 
tine -distance  in  the  two  layers.  The  tine-di  stances  in  the  upper  and 
lower  layers  ere  (OK  *  08 )ft1  end  dc/? 2>  reapeetivuly.  It  is  evident  from 
Figure  9  that  Od  -  CS  »  h/cos  a^j  also,  MS  •  x  -  2h  tan  a^. 

Bonce,  the  tine  for  the  refracted  wave  to  travel  path  OdCS  is 


Equation  6  any  be  simplified  by  replacing  the  trigonometric  functions  by 
their  equivalents  in  terms  of  the  velocities.  It  follows  from  Equation  5 
that 
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tan  a.  * 


Hence, 


■  *  V?T' 


The  t  las -distance  curve  far  the  refracted  wave  is  therefore  a  straight 
line  haring  a  elope  of  aagnitude  l/Vg  (figure  10). 
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The  tiaa-dietance  curve  can  also  be  used  as  a  basis  far  calculating  the 
depth*  h?  to  a  layer  of  a  higher  velocity.  Far  the  two  layer  case 
considered  Shove  It  can  be  found  that 
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The  critical  distance  x  is  the  distance  to  the  intersection  of  the 

c 

curves  for  the  direct  nave  and  the  refracted  wave  (beyond  the  range  of 
Figure  10).  This  expression  also  indicates  the  horizontal  distance  from 
the  source  to  t*u»  point  where  the  first  arrivals  from  the  higher  velocity 
layer  first  reach  the  ground  surface. 

la  the  usual  case*  the  geologic  section  consists  of  aore  than  two 
borlsootal  layers  of  thicknesses  by  ty  hy  ...  hQ  of  successively 
increasing  velocities  V^*  Vg*  V^>  ...  Vft  as  shown  in  Figure  11.  The 
paths  or  rays  panstratlng  to  the  various  layers  will  be  segnsnts  of 
straight  lines*  refracted  at  the  Interfaces  according  to  Snell's  lav.  The 
angle  of  the  path  taken  by  the  nlnlnaa  tins  ray  to  any  interface  can  be 
found  as  follows  (referring  to  Figure  11): 


AIWUT*-45-b7 

9nii  the  angle  Of  the  minimum  tir*  path  In  any  layer  is  determined  by  the 
velocity  in  that  layer  and  in  the  fastest  layer  penetrated  and  is 
indepeademt  of  the  velocity  in  intemsdiate  layers. 


Using  the  minimal  time  paths  given  above,  the  equation  for  the  time- 
distance  carve  far  the  refracted  oaves  can  be  expressed  in  the  following 
form: 


t  »  — *  + 
“  Vn 
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The  intercept  times  shown  on  the  top  of  Figure  11  can  be  written  as 
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Whenever  coincident-tine  curves  intersect  the  ground  surface ,  there  is  an 
abrupt  cnaage  in  slope  in  the  time --distance  curve,  indicating  that  at  such 
points  the  first  arrival  ia  due  to  a  different  set  of  waves.  As  the 
elopes  of  curves  of  a  single  system  are  not  discontinuous,  the  presence  of 
such  breaks  in  the  curve  is  a  signal  that  a  change  is  taking  place  from 
one  set  of  waves  to  another.  Unfortunately,  in  practice  it  is  not  always 
possible  to  secure  data  sufficiently  precise  to  detect  such  a  break. 


(3)  Partition  of  Energy  in  Reflected  sad  Refracted  Waves.  In 


tbs  previous  section,  the  lending  congressional  wavefront  was  traced  from 
its  source  throughout  the  surrounding  msdla,  including  adjacent  layers. 

It  was  sectioned  that  in  a  similar  esaner  the  first  arrivals  of  waves 
traveling  longer  paths  could  also  be  plotted,  although  the  problem  is 
complicated  if  the  seismic  velocity  has  been  altered  by  the  passage  of 
earlier  waves.  It  was  also  noted  that  the  isplngemsnt  of  either  &  con¬ 
gressional  or  distort! anal  wave  an  en  interface  separating  wlia  of 
different  characteristics  resulted  in  the  generation  of  four  new  waves. 
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To  gain  some  idea  of  the  nature  of  the  ground  action  at  a  point,  the 
contribution  of  all  naves  reaching  that  point  rust  be  evaluated.  The 
naves  reaching  the  point  by  direct  passage  through  a  single  medium  are 
usually  assumed  to  decrease  in  strength  with  distance  due  to  spatial  dis¬ 
persion.  At  the  interfaces  between  layers,  however,  it  is  necessary  to 
cospare  the  distribution  of  energy  among  the  reflections  and  refractions. 

Equation*  describing  the  partition  of  energy  between  the  reflected  and 
refracted  naves  at  an  Interface  separating  two  media  have  been  published 
by  Knott  (4)  and  are  available  in  more  recent  literature  (5).  numerical 
solutions  of  the  equations  for  a  vide  range  of  density  end  velocity  ratios 
were  published  by  Muskat  and  Msres  (6)j  by  Ihring,  Jsrdetaky,  and  Press  (7); 
also  by  Howell  (8)*  Because  of  their  application  to  this  study,  some  of 
the  results  of  the  work  by  Miskat  and  Meres  is  reproduced  here.  Selected 
results  from  their  companion  paper  (9)  are  also  presented. 

(A)  Energy  Division.  When  a  longitudinal  nave  impinges  on 
a  boundary  separating  layer*  of  different  elastic  constants  or  densities, 
its  energy  will  be  distributed  emcng  four  new  waves,  a  longitudinal  and  a 
transverse  wave  in  each  layer.  The  division  cf  energy  among  the  different 
waves  depends  on  the  velocities  and  densities  of  tbs  two  layers  and  on  the 
azgle  between  the  incident  wave  and  the  plans  of  the  Interface. 

To  illustrate  the  relative  importance  of  the  four  waves,  the  division  among 
0  them  of  the  energy  in  the  plane,  longitudinal,  incident  wave  is  tabulated 

fear  one  profile  in  Table  1.  The  ratio  of  the  density  of  the  medium  in  the 
lower  layer  to  that  in  the  upper  layer  is  1.1  and  Poisson's  rrtio  is 
assumed  to  be  0.29  for  both  madia.  The  parameter  *a”  is  the  ratio  of  the 
seismic  velocity  in  the  lower  layer  to  that  of  the  upper  layer. 

The  fractions  of  the  incident  longitudinal  wave  energy  converted  into 
energy  in  the  reflected  and  transmitted  longitudinal  and  transverse  waves 
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are  tiien  given  as  functions  of  ae isxLc  velocity  ratio*  a*  and  angle  of 
incidence  (vith  the  vertical)  9*  where: 

^  *  longitudinal  wave  reflection  coefficient 

T  -  transverse  wave  reflection  coefficient 
r 

*>  longitudinal  wave  transmission  coefficient 
-  Transverse  wave  transmission  coefficient 

The  sum  ef  the  coefficients  mat  equal  unity,  i.e., 

\  *  \  *  %  *  \  "  1 

It  nay  be  seen  from  Table  1  'that  at  angles  of  incidence  lees  than  critical 
a  large  fraction  of  the  energy  in  the  lac  idea?  wave  is  transferred  to  the 
refracted  longitudinal  wave.  Mach  of  this  transferred  energy  is  then 
refracted  into  still  lower  layers.  In  each  of  the  layers  into  which  waves 
are  refracted,  and  at  ranges  beyond  their  respective  critical  angles,  the 
refract  wave  propagates  outward  iu  a  direction  parallel  to  the  upper 
Interne-  c>i  %z*i  layer.  Uhder  these  conditions,  however,  the  transfer  of 
energy  across  the  Interface  1s  not  a  strong  one  and  only  »  nail  percen¬ 
tage  af  the  energy  in  the  refracted  waves  returns  to  the  upper  layer. 

It  is  evident  that  at  angles  of  incidence  greats*-  than  the  critical  value, 
all  of  the  energy  la  the  incident  longitudinal  wave  mat  be  reflected 
since  refractions  cannot  occur,  fbe  energy  flux  density,  or  energy  per 
unit  area,  in  the  reflected  wave  decreases  slightly  with  increasing  angles 
of  incidence,  reaches  a  minimum  and  then  Increases  slowly.  However  larger 
angles  of  incidence  Ug>ly  greater  ranges  and  for  cases  where  the  incident 
wave  is  spherical,  rather  than  plane,  spatial  dispersion  tends  to  offset 
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the  effect  of  the  increases  in  the  reflection  coefficients.  Nonetheless, 
the  net  result  is  that  the  reflections  have  a  greater  influence  on  the 
motion  in  the  upper,  softer  layer  than  the  refractions  returning  from  the 
lower  layer.  It  jwy  te  noted,  however,  from  Figure  8  that  the  first  wave 
arriving  at  the  surface  will  always  he  a  direct  wave  or  a  refraction, 
never  a  reflection. 

This  Conclusion  is  supported  by  the  work  of  Wolf  (10)  who  attempted  to 
explain  the  low  anplitudes  and  low  frequencies  of  refracted  waves  observed 
in  seismic  explorations.  Wolf  considered  a  profile  consisting  of  two 
semi -infinite  elastic  media  witn  both  the  source  of  disturbance  and  the 
recording  station  located  in  the  upper,  softer  medium.  In  a  numerical 
example,  he  assumed  that  the  seismic  velocity  in  the  upper  layer  was  10,000 
feet  per  second,  the  seismic  velocity  in  the  lower  layer  was  15,000  feet 
per  second?  and  the  range  was  5000  feet.  The  disturbance  was  a  displace¬ 
ment  pulse  of  0.02  second  duration. 

The  calculated  horizontal  displacement  at  the  recording  station  is  shown 
in  Figure  12.  Kie  refracted  wave,  RR,  arrived  at  the  station  at  time  0.33 
seconds  after  l he  initiation  of  the  disturbance  and  consisted  of  a  unidi¬ 
rectional  pulse  which  lasted  until  time  0.5  seconds  when  the  direct  wave, 

P,  arrived.  Shortly  after  the  arrival  of  the  direct  wave,  the  motion  due 
to  a  reflection,  RF,  from  the  interface  was  observed. 

Wolf  notes  that  for  this  case  the  amplitude  of  the  horizontal  component  of 
the  displacement  due  to  the  refrectio  wee  only  0.08  timej  the  amplitude  of  the 
direct  wave  and  that  the  direction  of  the  returning  refracted  wavefront 
was  that  of  total  reflection,  i.e.,  at  the  angle  of  critical  refraction. 
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Figure  12.  Waveform  of  Horizontal  Ground  Motion: 

RE  -  Refracted  Wave*  D  -  Direct  Wave* 
RF  -  Reflected  Wave  (10). 


The  previous  discussion  has  been  concerned  with  waves  generated  by  a 
single  point  source  of  energy.  However  the  nuclear  burst  also  generates 
a  strong  airblast  wave  which  loads  the  ground  surface  progressively  as  it 
propagates  outward  from  the  origin.  As  a  result*  the  strengths  of  the 
refracted  waves  could  be  reinforced  by  the  loading  of  subsequent  points 
as  the  air  wave  advances  over  the  ground  surface.  The  direct  and  reflec¬ 
ted  waves  could  also  be  reinforced  in  the  sane  manner.  The  degree  of 
reinforcement  would  also  depend  upon  the  seismic  velocities*  thicknesses 
of  the  layers*  and  velocity  of  the  air  wavefront.  Since  both  the  reflec¬ 
ted  and  refracted  waves  result  from  the  incidence  of  the  direct  wave  on 
an  interface*  it  is  probable  that  if  reinforcement  were  to  occur  from 
the  air  wave*  all  waves  would  benefit  nearly  proportionately  and  the 
amplitude  ratios  would  remain  about  the  same.  This  reasoning  will  be 
examined  further  in  the  section  comparing  the  actual  test  data  with  the 
predictions  using  the  wavefront  diagram. 

The  transmission  and  reflection  coefficients  contained  in  Table  1  are 
derived  for  a  two  layer  case.  The  results  may*  of  course*  be  applied  to 
a  multiple  layer  system  in  a  step-by-step  proces  .  However,  in  order  to 
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•hew  the  influence  of  the  eagle  of  incidence  and  the  Influence  of  multiple 
•trail  fleet  ion,  Musfcat  and  Meres  (9)  presented  a  number  of  examples, 
figure  13  is  taken  frost  their  work. 

gxaadnfttlon  of  the  examples  of  Figure  13  show  that,  when  the  angle  of 
incidence  is  less  than  the  critical  angle,  the  effect  of  the  magnitude  of 
the  angle  of  incidence  is  small.  For  most  estimates  of  the  energy  return, 
an  assuaptlon  of  normal  incidence  is  Justified,  and  the  angle  of  incidence 
with  the  deepest  layer  should  suffice  for  all  purposes  when  normal  Inci¬ 
dence  is  not  explicitly  assumed. 

It  can  also  he  seen  from  the  examples  that  the  presence  of  intermediate 
interfaces  greatly  attenuates  the  energy  ultimately  returning  from  the 
layer  of  greatest  depth.  A  single  intermediate  interface  reducing  the 
total  velocity  contrast  by  two  will  reduce  the  fraction  of  incident  energy 
returning  from  the  highest  velocity  layer  by  a  factor  of  five.  An  addi¬ 
tional  Intermediate  Interface  causes  a  further  reduction  by  a  factor  of 
about  five,  so  that  the  energy  returning  ultimately  from  the  16,000  feet 
per  second  velocity  layer  would  be  only  0.4  percent  in  Figure  13(g),  as 
compared  to  11.1  percent  of  the  incident  energy  in  Figure  13(a).  Part  of 
this  attenuation  is  due  to  the  energy  reflected  from  the  intermediate 
interfaces.  However,  about  half  of  the  total  decrease  is  due  to  the 
higher  ultimate  loss  of  the  incident  energy  by  transmission  into  the 
highest  velocity  layer  because  of  the  reduced  velocity  contrast. 

Figures  13(h)  and  I3(i),  uhich  correspond  in  their  gross  features  to  the 
type  of  section  encountered  in  certain  Midwestern  districts,  show  again 
the  high  degree  of  attenuation  resulting  from  the  intermediate  reflection 
processes.  Under  the  assumed  conditions,  less  than  one  percent  of  the 
incident  energy  will  return  from  the  highest  velocity  and  deepest  layer, 
whereas  more  than  five  percent  will  return  from  the  two  upper  interfaces. 
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The  velocity  section*  of  Figure*  13U)  and  13 (k)  show  the  influence  of  an 
intermediate  lev -velocity  layer.  It  is  seen  that,  if  the  section  resumes 
its  normal  velocity  increase  with  depth  after  passing  the  icw-velocity 
layer,  the  high  contrast  at  the  base  of  the  latter  will  give  rise  to  high 
reflection  coefficients  and  greater  ultimate  returns  of  reflected  energy  to 
the  surface  than  viU  Motions  where  the  velocities  increase  monotonically 
wlth  depth  and  where  all  the  contrasts  are  rather  small.  Thus,  in  the 
specific  cases  shown  in  Figures  13(d)  and  13 (k),  almost  as  much  energy  will 
be  returned  from  the  deepest  layer  as  from  the  first  interface.  Moreover, 
in  coopering  Figure  l3(j)  with  13(h),  more  than  three  times  as  much  energy 
will  be  reflected  from  the  15 >000  feet  per  second  layer  in  Figure  13  (j ) 
than  in  Figure  13(h),  although  the  latter  has  one  less  intermediate 
interface. 

The  percentage  of  the  energy  reflected  from  a  weathered  or  indistinct 
interface  between  layers  would  be  significantly  reduced  below  the  values 
given  in  the  tables  or  examples.  An  indistinct  boundary  between  layers 
of  a  geologic  formation  is  often  the  result  of  a  gradual  change  in  the 
conditions  of  deposition  or  weathering  prior  to  deposition  of  the  suc¬ 
ceeding  layer.  In  this  case,  the  velocity  gradient  is  gradual  and  may 
result  in  more  of  the  energy  being  refracted  into  the  deeper  formation 
and  less  being  reflected. 

The  percentage  of  energy  reflected  at  a  discontinuity  Increases  rapidly 
as  the  angle  of  incidence  approaches  the  critical  refraction  angle. 

However,  Gutenberg  (ll)  showed  that  there  would  be  no  corresponding 
Increase  in  the  amplitudes  of  the  reflected  waves  at  the  ground  surface. 
Instead,  the  maximum  amplitudes  of  reflected  waves  are  found  near  the 
source  of  the  disturbance.  This  is  because  the  ground  motion  is  due  both 
to  the  reflection  from  the  interface  and  to  its  subsequent  re-rtflection 
downward  from  the  ground  surface.  The  amplitude  also  depends  on  the  rate 
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of  change  of  the  angle  of  incidence  with  distance.  That  is,  the  larger 
the  angle  of  incidence  the  larger  the  area  will  he  at  the  ground  surface 
over  which  the  initial  energy  Intensity  will  be  spread. 

(4)  Use  of  Time -Distance  Curves  and  Wavefront  Diagram  with 

Airblast  Arrival  Curve.  The  time-distance  curves,  whether 

derived  from  the  construction  of  a  wavefront  diagram  or  plotted  from  a 

seismic  survey,  can  be  used  in  conjunction  with  the  airblast  arrival  curve 

to  predict  the  point  of  first  outrunning  and  the  arrival  time  of  various 

ground  waves  at  the  point  of  Interest.  Figures  lb  through  17  present  the 

time-distance  curves  derived  from  the  wavefront  diagram  of  Figure  5,  and 

the  airblast  arrival  curve  for  a  20  Ml  weapon,  (Figure  16).  Hie  series  of 

figures  illustrates  the  procedure  used  in  locating  the  range  where  the 

arrival  time  of  the  refracted  wave  coincides  with  that  due  to  the  airblast. 

At  greater  ranges  the  refracted  wave  will  arrive  prior  to  the  airblast 

pressure.  In  addition  to  the  arrival  times  of  the  refracted  waves,  the 

arrival  times  of  the  direct  wave  and  the  reflected  wave  from  the  two 

deeper  layers  are  also  indicated. 

Hie  range  at  which  outrunning  first  occurs  nay  be  found  by  eliding  the 
origin  of  the  refracted  wave  time-distance  curve  along  the  airblast 
arrival  time  curve  until  the  two  Intersect  at  the  minimum  ground  range. 

This  point  will  be  found  when  the  origin  of  the  refraction  time-distance 
curve  is  located  at  the  point  cm  the  airblast  arrival  curve  where  the 
velocity  of  the  airblast  wave  equals  the  velocity  of  the  layer  whose 
time-distance  curve  intersects  the  airblast  wave  arrival  time  curve.  For 
example,  in  Figure  lb  the  origin  is  located  at  a  point  on  the  airblast 
curve  where  the  airblast  front  is  traveling  with  a  velocity  of  approxi¬ 
mately  10,000  feet  per  second.  Ihls  is  equal  to  the  velocity  of  layer  3* 

In  this  example,  outrunning  occurs  first  on  the  ground  surface  at 
Station  6200.  Hie  origin  of  the  refracted  wave  arriving  first  at 
Station  6200  was  the  disturbance  generated  by  the  airblast  wave  when  it 
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Figure  16.  Location  of  Range  of  First  Outrunning  with  Source  of  a  Range  of  5000  Feet 
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was  at  a  range  of  3300  feet.  The  peak  overpressure  at  Station  6200  is 
300  pounds  per  square  inch*  and  at  Station  3300  is  1300  pounds  per  square 
inch  (27), 

It  any  be  noted  in  Figures  14  through  17  that  the  refracted  wave  reaches 
stations  betveen  6200  and  6300  feet  when  the  origin  varies  over  a  range 
uetveen  2500  and  3000  feet.  Due  to  the  velocity  variation  of  the  airblx^t 
wave  and  the  seismic  velocities?  all  the  energy  carried  by  the  critically 
refracted  waves  between  this  range  will  arrive  at  nearly  the  sane  tine. 

Figures  t Jsrn^i  23  indicate  the  contours  on  which  disturbances  gener¬ 

ated  m  tiiSK-  t  arrive  simultaneously  at  the  point  of  interest.  Figure  19 
is  a  plan  view  of  the  ground  surface  and  contains  one  set  of  concentric 
curves  radiating  from  Ground  Zero  (GZ)>  representing  the  position  of  the 
air  blast  wavefront  at  various  tine  intervals.  The  airblast  data  are  fron 
Operation  Tumbler  Shot  Ho.  1,  and  are  presented  with  the  wavefront  diagraa 
in  Figures  24  and  25.  Tne  concentric  circles  radiating  fron  Station  5V 
are  lines  of  constant  transit  tins  fron  a  disturbance  on  the  surface  to 
Station  5V.  The  transit  tines  were  obtained  fron  a  modified  tine-distance 
curve  (Figure  26)  for  Frenchman  Flat  (FF- 1 ) .  The  tine-distance  curve 
for  this  site  in  turn  was  obtained  by  construction  of  a  refraction  wave- 
front  diagram  (Figure  27).  The  somewhat  elliptical  contour®  were  then 
drawn  through  points  which  yielded  equal  transit  times  from  Ground  Zero 
to  Station  5V.  The  diagrams  and  figures  are  used  later  to  coapare  the 
predictions  of  the  direction  of  motion  and  arrival  tines  with  the  field 
gp  data  fro*  the  shots  performed  at  this  site  durir$  Operation  Tumbler  Shot 

Ho.  1. 

It  can  be  notad  by  examination  of  Figure  19  that  the  areas  within  which  an 
impetus  of  finite  duration  may  be  generated  end  My  arrive  at  the  station 
at  the  same  time  represent  a  wide  range  of  instantaneous  pressures  awl 
transit  distances.  In  the  particular  condition  used  to  construct 
Figure  19*  it  can  be  seer.,  for  enable ,  that  the  refracted  wave  at  an 
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arrival  time  of  0  Vj  seconds  is  generated  from  both  the  far  side  of  Ground 
Zero  and  from  points  beyond  tne  station  of  interest.  Tht  same  general 
situation  is  true  for  the  arrivals  of  the  other  waves  as  shown  on  the 
other  fjgures.  This  is  one  of  the  several  conplications  which  make  the 
prediction  of  energies  at  a  given  station  so  difficult  by  graphical 
methods.  General  observations  regarding  the  intensities  of  the  various 
arrivals  will  be  discussed  in  the  section  conparing  the  predictions  with 
the  actual  field  data. 

Figure  2k  represents  the  position  of  the  refracted  wavefront  when  the 
refracted  wavefront  diagram  shown  in  Figure  2 7  is  combined  with  the 
airblast  arrival  data  given  in  Figure  25-  It  differs  from  the  refraction 
diagram  of  Figure  27  in  that  in  Figure  27  the  refracted  wavefronts  emanate 
from  a  point  source  at  "0"  while  in  Figure  25;  the  refracted  wavefronts 
are  generated  by  the  airblast  wave  moving  alonj  the  ground  surface. 

Figure  25  is  constructed  by  shifting  the  origii  of  the  wavefront  diagram 
(Figure  27)  so  that  it  represents  the  position  of  the  airblast  front. 

The  various  positions  of  the  refracted  wave  are  sketched  and  their  tran¬ 
sit  time  recorded.  The  envelopes  which  represent  the  most  advanced 
positions  of  the  wavefront  at  a  given  time  are  those  shown  cm  Figure  2k. 
The  times  shown  represent  the  total  elapsed  time  since  the  burst.  Since 
Operation  Tuafeler*  Shot  No.  1  was  an  airblast.  the  initial  time  at  Ground 
Zero  wao  equal  to  0.52  seconds. 

Figure  2k  shews  the  general  change  of  direction  and  slope  of  the  first 
arrival  wavefronts  as  they  recede  from  the  point  of  disturbance.  From 
this  curve  it  Is  possible  to  scale  the  arrival  times  and  direction  of 
the  first  refracted  wave  at  any  point  below  the  surface .  This  diagram 
does  not  show  the  area  from  which  a  given  wave  was  generated.  If  desir¬ 
able,  it  would  be  possible  to  construct  in  a  similar  amoner  a  diagram 
showing  the  propagation  of  reflected  waves. 
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C.  Comparison  of  Predictions  vith  Field  Data.  In  this  section  the 
acceleratio^-tia*  waveforms  of  several  test  shots  are  compared,  with 
predictions  based  on  the  wavefront  diagram  method  outlined  in  this  report. 
Field  conditions  and  measurements  used  in  this  comparison  have  been 
obtained  from  the  unclassified  reports  of  Operation  Tuablea,  Shots  Ifo.  1 
through  No.  4  (12)  and  Iron  Operation  Flumbbob,  Shot  Priscilla  (15).  Data 
from  Operation  Hardtack- 1  (ik)  are  not  presented  due  to  their  classified 
nature;  howeve . ,  &<..  ?  of  the  pertinent  observations  obtained  from  a  com- 
parish  a  w  .th  the  date  from  those  shots  will  be  discussed. 

Pirst>  seismic  velocities  computed  from  data  obtained  in  a  refraction 
survey  at  the  test  site  are  compared  with  those  indicated  by  the  shoe': 
arrival  times  during  a  nuclear  test.  The  comparison  of  velocities 
Measured  during  a  low  stress  seismic  survey  vith  those  of  the  much  higher 
etreos  test  should  indicate  the  variation  of  ways  velocity  with  wave 
strength.  A  considerable  volume  of  material  is  available  concerning  the 
theoretical  and  experimental  inelastic  and  dynamic  behavior  of  soils  (i5)» 
therefore,  a  discussion  of  this  aspect  will  not  be  included  in  this  paper. 

In  the  second  part  of  this  section  a  comparison  Is  made  of  the  predicted 
arrival  times  and  directions  of  motion  using  the  wavefront  diagram  method 

with  the  field  test  records.  ..w  ’  * 

..  y.  .  y  ' 

(l)  Field  Test  Data.  Field  test  data  used  in  tha  remaining 
portions  of  this  section  are  presented  in  Bsbies  2  througn  12. 
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Table  2 

miS  TEST  BATA,  OPERATION  TUMBLES  (12) 


LOCATION 


1  April  1952 
15  April  1552 
22  April  1Q52 
1  May  1952 


HEIGHT 
OP  BURST 
(ft) 


YIEIS 

(KT) 


1.05 

1.15 

50.00 

19.60 


Table  5 

afittftUC  VELOCITIES  (MODIFIED);  OPERATION  TUMBLER  SIZES  (12) 


VKaCHMr  FIA* 

YUCCA  FIAT 

T-7  area 

— 

SEISMIC 

SEISMIC 

VELOCITY 

VELOCITY 

(ft) 

(fP») 

<">  _ 

(fpe) 

\ 

0-20 

1500 

0-5 

800 

20-60 

2A00 

1  5-100 

VARIABIE 

(  .200  fps  5 BED) 

60-200 

2800 

ICC *2, 5 

3500 

200-600 

5600 

275-300 

5COO 

600 

10,300 

500 

600C 

OPERATION  TUMBUER ,  SHOT  NO.  1  -  GROUNDSHOCK  DATA 
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Table  10  ■ 

''PERATI0H  TUMBLER,  SHOT  10.  3  -lAIFBl'iST  DATA 
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p  (bflt  IB  has  three  separate  shocks; 

arrival  tlaes,  0.28V,  0.29$,  and  0.307  sec; 
pressures,  26,  V2,  sad  V$  psi. 

q  Record  not  clear;  possibly  a  precursor  exists  vith  arrival 
tins  of  0.2V$  sec.,  aad  pressure  ?+  <.2  pal. 
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(2)  Comparison  of  Seismic  Velocity  vlth  Test  Velocity.  A*  a 
wfiu  of  compering  the  velocity  measured  luring  a  normal  refraction 
sciatic  aur/ey  vlth  the  velocity  of  the  groundshock  v-.ve  from  a  nuclear 
explosion,  time -distance  curves  were  constructed  using  the  published 
seismic  velocities  of  the  test  sites.  The  time -distance  curve  vas  then 
shifted  along  the  aixblast  arrival  curve  until  the  point  of  outruztJng  vas 
determined,  the  results  of  this  operation  are  presented  in  Figures  29 
through  JJ. 

9m  coaperisan  of  the  predicted  times  of  arrival  of  the  outrunning  ground 
eaves  with  thoee  actually  observed  for  the  Tumbler  series  of  shots  shews 
rssstrkajole  agraewnt.  It  should  be  noted  that  the  range  of  overpressures 
for  the  TttSbler  series  Is  relatively  low,  varying  between  136  and  2  pounds 
per  square  inch-  T3iu* ,  at  these  pressures,  the  seismic  survey  velocities 
appear  to  be  good  approximations  to  use  in  predicting  first  arrivals  for 
outrunning. 

9m  comparison  shown  in  Figure  90,  using  the  Shot  Priscilla  data,  was  not 
conclusive  since  outrunning  did  not  occur  at  the  test  stations.  A  com¬ 
parison  of  the  velocities  calculated  from  the  arrival  of  the  peek  pressure 
nave  at  the  gage*  between  tbs  eurfsee  and  a  depth  of  90  feet  at  Frenchmen 
Fist  at  yresvrnbly  the  100  pentads  per  square  lash  level  and  thoee  measured 
*  refraction  seismic  survey  la  reported  by  Sauer  (U ).  9m  velocity  of 
the  petit  stress  save  is  shown  to  be  between  2/9  sad  3A  of  that  measured 
by  the  reiemdc  survey. 


A  survey  aede  after  Shot  Priscilla  showed  no  permanent  dioplacesmnt  below 
nvarprseaurss  of  about  ho  pounds  par  square  inch.  Persmnent  displacement 
of  coarse,  would  imply  Inelastic  behavior  of  the  soils. 

Id  contrast  to  the  dry,  leesrt  type  plays  soils  of  Frenchman  and  Yucca 
Flats  at  the  Pseuds  Vast  Site,  the  soil  conditions  of  the  Pacific  Proving 
Oronede  are  saturated  coral  anode  sod  coral  conglnaerata.  Epical  of  the 
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coral  atoll  of  the  Biivetok  Proving  Qrom»ds,  site  of  Shot  Koe  and  Cactus, 
Is  a  surface  layer  5  to  7  feat  thick  of  dry,  loose  coral  sand  and  broken 
■hello.  The  seismic  velocity  of  the  surface  layer  is  approximately  800 
feet  per  second.  The  surface  layer  is  underlain  by  the  water  surface, 
and,  in  places,  by  a  very  bard  layer,  several  feet  thick,  of  cemented 
coral  conglomerate.  The  seismic  velocity  of  the  conglomerate  is  on  the 
order  of  8,000  feet  per  second.  Belov  the  water  table  the  sites  are 
underlain  by  alternating,  random,  and  discontinuous  layers  of  loose  to 
dense  coral  sand  and  shells,  and  other  remnants  of  cemented  coral  con¬ 
glomerate  shells.  The  velocity  of  the  saturated  sands  is  on  the  order  of 
5,000  feet  per  second. 


t 


*»»»*» 


Figpars  3?.  Oaapariaoa  att  Yield  Data  with  Calami Mod  tiam-Dlatane* 
Curve,  Operation  fsabler.  Shot  Wo.  l. 
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Figure  32.  Comparison  of  Field  Date  with  Calculated  Tirae-Di  stance 
Curve,  Operation  'Ambler,  Shot  No.  3. 


Figure  33.  Coapariaon  of  Field  Data  with  Calculated  ttme-Di stance 
Curve,  Operation  flakier,  Shot  No.  4. 
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The  test  shot  for  Koa  was  arranged  so  that  the  first  gages  recording 
groundshock  were  located  at  stations  where  the  overpressures  were  about 
1000  pounds  per  square  inch.  Gages  were  also  located  at  stations  which 
would  receive  overpressures  of  about  100  pounds  per  square  inch. 

®e  significant  feature  of  the  ground  motions  for  Shots  Koa  and  Cactus  was 
that  the  waveforms  had  the  character  of  the  refracted  wave,  but  of  large 
amplitude,  (Figures  3  and  k).  The  velocities  of  these  very  strong  waves 
were  slightly  higher  than  measured  by  the  refraction  seismic  surveys  at 
the  respective  sites  (l4). 

It  is  probable  that  the  hard  coral  conglomerate  acted  nearly  elastically, 
as  perhaps  did  the  saturated  coral  sands.  Although  the  final  result  of 
the  pressure  and  vibrations  on  the  saturated,  lc"se  coral  sand  was  deusi- 
fication,  this  may  not  have  occurred  until  the  water  had  time  to  drain 
from  between  the  sard  grains.  As  a  result,  the  sand-water  system  may 
have  acted  elastically  during  the  passage  of  the  first  high  intensity 
waves.  This  possibility  may  also  explain  why  the  maximum  displacement  did 
not  occur  until  a  time  after  the  entire  airwave  had  passed  over  the  recording 
stations. 

The  inelastic  behavior  of  a  soil  wvjuld  he  most  evident  when  e  large  change 
ia  stress  occurs  relative  to  the  initial  stress  of  the  soil.  Considering, 
for  a  typical  site  that  is  not  saturated,  that  the  vertical  stress  due 
the  overburden  is  about  0.9  pounds  per  square  inch  per  foot  of  depth, 
at  depths  on  the  order  of  100  to  200  feet,  the  soil  is  initially  stressed 

to  and  under  equilibrium  at  pressures  of  90  to  180  pounds  per  square 
inch.  The  overpressures  at  the  ground  surface,  considering  attenuation  of 
■tress  with  depth,  may  havs  to  be  such  higher  than  the  initial  stress  of 
the  soil  before  the  velocity  of  the  ground  waves  is  significantly  reduced 
by  nonlinearities. 
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The  velocity  at  which  soils  below  the  water  level  transmit  a  high  stress 
wave  would  be  close  to  that  of  the  seismic  velocity  of  water.  Thus  the 
velocity  of  a  saturated  soil  as  measured  during  a  seismic  survey  would 
seem  to  be  a  good  approximation  to  use  in  constructing  a  wavefront 
diagram. 


Sauer  (16)  suggests  that  in  predicting  peak  amplitudes  of  ground  motion 
three  fourths  of  the  seismic  velocity  for  soils  and  fractured  rock  should 
be  used.  Based  upon  the  evidence  shown  by  the  comparison  of  the  predicted 
arrival  time  of  the  refracted  wave  and  that  actually  measured  during  the 
weapon  tests ,  it  would  seem  reasonable  to  use  the  three  fourths  reduction 
for  dry  soils  at  overpressures  over  40  pounds  per  square  inch;  and  for  the 
case  of  saturated  soils ,  rock,  and  perhaps  soil  at  depths  great--"  than  200 
feet ,  to  use  the  velocity  as  measured  by  a  seismic  survey. 

(3)  Comparison  of  Predicted  Arrival  Times  and  Direction  of 
Moiion's  with  Test  Acceleration  Waveforms .  To  demonstrate 

the  manner  in  which  the  contours  of  equal  arrival  time  (Figures  19  through 

23)  may  be  combined  with  the  time -distance  curves  (Figures  26  through  28} 

to  yield  arrival  times  and  directions  of  motion ,  and  aa  an  indication  of 

the  validity  of  the  predictions*  test  data  fron  Operation  Tumbler,  Shot 

No.  1,  are  compared  with  the  predicted  motions  in  Figures  34,  33  and  36. 

At  the  top  of  each  figure,  acceleration  waveforms,  recorded  at  Stations  3 

and  5>  are  reproduced  here  from  Reference  12.  The  station  number  and 

direction  of  motion  are  indicated  by  the  symbols  near  the  vertical  scale; 

that  is,  "3V"  denotes  that  the  record  was  taken  at  Station  3  aid  that  the 

vertical  motion  was  measured.  The  letter  "H"  indicates  horisontal  motion 

while  tangential  motion  is  labeled  "T".  the  lower  charts,  characterised 
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by  the  arrows,  show  the  predicted  arrival  tins  of  the  refracted,  direct 
and  reflected  waves  and  their  directions  of  motion. 

flnhesfjnnnt  to  tLe  earliest  arrival  of  each  wave  indicated  by  the  arrows, 
tbe  cotlon  continues  as  long  as  the  disturbance  reoalns.  The  position  of 
the  arrows  after  the  first  arrival  therefore  is  somewhat  arbitrary  and  is 
indicative  only  of  the  general  notion  at  these  later  times,  lbe  intensity 
of  the  notion  varies  as  the  intensity  of  the  source  changes.  However,  the 
lengths  of  ths  arrows  shown  in  the  figures  are  not  intended  to  indicate 
relative  strengths  of  tbe  waves. 


The  time  of  arrival  of  the  first  of  each  type  of  wave  can  be  determined 
directly  fro*  the  contours  of  figures  19  through  23.  Knowing  the  range 
from  the  source,  the  direction  of  notion  of  refracted  waves  will  be  normal 
to  the  refraction  wavefront  at  that  range  as  shown  in  Figure  27.  From 
Figure  28,  the  direction  of  notion  of  reflected  waves  can  be  found  in  a 
similar  manner. 


If  contours  of  equal  arrival  tine  have  not  besot  plotted  for  the  site  and 
weapon  conditions,  tbs  sans  results  can  be  obtained  fron  time-distance 
curves  such  as  that  shewn  in  Figure  26.  This  figure  was  exqplpyed  In 
predicting  the. Arrival  times  shown  in  Figures  3^  and  35  far  Station  3. 

sliding  the  origin  of  the  time-distance  curve  along  the  airblast 
arrival  curve  in  the  same  manner  that  tbe  point  of  outrunning  is  located, 
the  total  time  required  for  each  wave  to  travel  a  given  range  nay  be 
determined.  It  will  be  found  to  be  faster,  however ,  to  draw  contour  naps 
of  equal  arrival  time.  In  either  case,  tbs  direction  of  motion  Is  ob¬ 
tains*  by  the  use  of  the  appropriate  wavefront  diagram  as  before.  It 
should  he  noted  that  time-distance  curves  must  be  constructed  for 
horisoctal  planes  passing  through  tbe  point  of  interest. 

The  geologic  model  used  in  this  cotqparlsoc  was  modified  somewhat  from  the 
test  data  in  order  to  simplify  the  calculations  and  in  an  attempt  to 
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correct  for  the  lower  effective  velocity  of  the  surface  layer  while  loaded 
by  the  overpressure.  The  velocity  profile  used  appears  on  Figure  28.  In 
reality ,  the  seismic  velocity  profile  measured  in  a  boring  by  Shannon  and 
Wilson  (3)  indicated  the  existence  of  a  layer  with  a  velocity  of  2500  feet 
per  second  sandwiched  between  the  two  surface  layers  with  a  velocity  of 
1500  feet  per  second.  As  shown  in  Figure  13>  when  such  a  condition  exists 
the  reflected  energy  from  the  next  deeper  high-velocity  layer  is  nearly  as 
high  as  that  from  the  first  high-velocity  layer. 

The  presence  of  this  2500  feet  per  second  layer  may  “''count  for  the 
apparently  late  reflection  arrivals,  since  this  layer  was  not  included  in 
the  construction  of  the  wavefront  diagram  and  other  curves. 

There  appears  to  be  reasonably  good  agreement  between  the  arrival  times 
and  between  the  phasing  of  the  vertical  motion  with  the  horitontal  motion. 

It  appears  quite  clear,  in  the  cases  shown,  that  the  first  arrivals  of  the 
refraction  wave  do  not  carry  much  energy  in  relation  to  the  direct  ~#ave 
aid  early  reflections.  At  later  times,  when  energy  is  arriving  from  many 
sources,  the  waveform  could  be  expected  to  have  a  widely  fluctuating 
pattern.  Relatively  minor  chaiges  in  the  sequence  of  stratification  would 
significantly  alter  the  pattern. 

The  negative  (downward)  acceleration  oc curing  near  t  -  0.8  seconds 
(Figure  36)  is  not  indicated  in  the  diagram  aa  being  initiated  by  the 
arrival  of  a  wave.  Hot  all  waves,  of  course,  have  been  included  in  the 
diagram.  Shear  waves,  plastic  shock  fronts,  and  reflections  due  to  the 
iqpiigemeat  of  the  refracted  waves  on  the  ground  surface,  for  example, 
have  not  been  represented  on  the  time -distance  curve  although  some 

greater  detail  could  have  been  shown  If  desired. 

* 

It  is  more  likely,  however,  that  the  fall  cycle  oscillation  of  the  accele¬ 
ration  curve  is  due  to  the  relatively  short  duration  of  the  airblast.- loading 
which  generated  the  refraction.  At  the  range  at  which  the  refracting  wave 
was  initiated,  the  time  of  raying  of  the  overpressure  to  fifty  percent  of 
its  peak  value  was  less  than  o.io  seconds.  The  rapid  decay  of  the  pressure 
would  propagate  as  a  family  of  expansions  which  would  account  for  the  oscil¬ 
lation  but  wt>i«h  have  not  been  shown  on  the  time  distance  curve. 
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Although  the  first  reflections  from  lower  layers  to  arrive  at  the  station 
cone  from  sources  varying  widely  in  range  and  overpressure,  the  work  of 
Gutenberg  (ll)  shows  that  the  close-in  sources  contribute  the  largest 
percentage  of  energy  even  though  they  nay  not  represent  reflections  of 
rays  flatter  than  the  critical  angle. 

Ho  quantitative  estimates  can  be  node  regarding  the  energy  associated  with 
the  various  waves.  The  evaluation  of  the  stress  at  a  point  by  spatial 
dispersion  using  the  relationships  of  Boussinesq  or  Vestergaard  have  been 
applied  to  given  areas  upon  which  the  load  intensity  varies  with  tine . 

This  approach  nay  be  feasible  for  gross  low-frequency  characteristics  of 
the  waves,  but  probably  would  not  reveal  the  detail  needed  in'. this  appli¬ 
cation. 

D.  Jonclualona.  The  wavefront  diagraa  is  a  useful  and  staple  aeons 
far  predicting  certain  wave  fora  parameters  of  importance  to  shock  isolation 
system  design.  If  the  site  geology  is  known  in  reasonable  detail,  the 
type  of  wave  first  arriving  at  the  point  of  interest,  the  first  arrival 
times  of  waves  due  tu  the  airblast,  those  transmitted  directly  from  the 
source,  and  those  reaching  the  point  of  interest  by  reflections  or  refrac¬ 
tions  from  umier lying  layers  can  be  estimated  with  good  accuracy.  Sms, 
the  time  differences  between  the  arrivals  of  the  various  waves  can  be 
determined. 

lfonllnearltles  in  the  stress-strain  relationship  of  the  soil  are  of  little 
importance  in  calculating  the  arrival  time  of  the  first  wave  to  reach  the 
point  of  interest,  as  the  wave  will  travel  at  the  seismic  velocity  of  the 
undisturbed  medium.  for  subsequent  waves,  it  is  found  that  the  initial 
state  of  stress  at  depths  where  high-intensity  waves  from  large  yield 
weapons  would  be  reflected  or  refracted  is  usually  sufficiently  high  so 
that  the  soil  modulus  remains  fairly  constant.  Therefore,  ...  ~egione 
outside  the  rupture  tone,  it  is  probable  that  the  waves  will  travel  at 
velocities  nearer  to  the  measured  seismic  velocity  than  to  those  which 
might  be  predicted  sn  the  basis  of  an  initially  unstressed  medium. 
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On  tbe  teals  of  ccsqjerisons  of  the  predicted  arrival  time*  of  refracted 
nm  irlth  those  Measured  during  several  weapons  tests*  it  appears  to  be 
reasonable  to  base  calculations  on  a  wave  velocity  equal  to  three  fourths 
of  tbs  Measured  selsalc  velocity  for  dry  aoils  and  pressures  in  excess  of 
MO  pounds  per  square  Inch.  For  saturated  soils  and  rock  and  for  all  types 
of  soil  forastlons  at  depths  greater  than  200  feet*  the  use  of  the 
mbs  mured  sslaaic  velocity  appears  to  yield  good  correlations  with  KS  and 
00  data. 

In  Reference  1,  it  vas  suggested  that  tbe  type  H  wave  was  tbe  effect  of 
an  Interaction  of  the  alrblast  wave  with  waves  refracted  from  lover 
surfaces.  However  the  present  study  has  shewn  that  relatively  snail 
amounts  of  energy  are  carried  back  into  the  upper  layers  by  refracted 
waves.  Instead*  it  Is  the  reflections  which  are  primarily  responsible 
for  the  return  of  energy  froa  s  lewer  to  an  upper  layer. 

for  this  reason  it  is  recommended  that  the  nornal  refraction  survey  nade 
during  site  evaluation  be  supplemented  with  a  reflection  seismic  survey 
of  the  arte  near  the  proposed  ft  :illty  location.  The  refraction  survey 
will  yield  data  which  can  be  used  directly  in  the  construction  of  the 
wavefront  diagram  as  described  earlier  in  this  section.  Date  from  the 
reflect ion ' selsbogreas  will  not  only  be  of  value  in  constructing  tbe 
wavefront  diagram  but  rtould  also  assist  in  evaluating  the  relative 
asgUtudes  of  the  reflections  from  lower  layers. 

farther*  while  the  alrblast- induced  wave  may  contribute  to  the  oscillatory 
nature  of  the  type  II  wave,  combinations  of  refractions*  and  reflections, 
with  the  direct  ccmpressicoal  and  expansion  waves  can  produce  oscillations 
without  the  necessity  fur  the  premier  of  the  alrblast- induced  component. 

from  the  wavefront  diagram  and  supporting  calculations  the  direction  of 
each  wave  reaching  the  point  of  Interest  from  each  of  tbe  sources  can  be 
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determined  and,  with  certain  assumptions,  the  energy  in  each  ray  grossly 
approximated.  It  is  evident  from  the  proceeding  discussions,  however, 
that  the  total  number  of  waves  emanating  from  eac.  source  is  very  large. 

In  the  exasqples  only  the  simplest  of  the  coapresslonal  body  waves  have 
been  considered,  the  loading  function  has  been  assumed  to  be  invariant 
with  time,  and  shear  and  surface  waves  have  been  neglected  entirely. 

While  propagation  of  all  of  these  phenomena  could  be  included  in  the 
time -distance  curve  and  wavefront  diagrams  could  be  drawn,  the  fact  that 
only  the  roughest  estimates  can  be  made  of  the  time-history  of  the  motion 
djue  to  each  of  the  waves  severely  restricts  the  usefulness  of  this 
graphical  method  for  predicting  the  waveform  of  the  Type  JJ  ground  motion. 

Nonetheless,  the  wavefront  diagram  and  the  time  d i stance  curves  can  be 
useful  tools  in  providing  «  simple  means  for  better  visualising  the 
routes  by  which  certain  waves  reach  the  point  of  interest  ard  for 
determining  quantitatively  sow  ijportant  time  parameters  of  the  ccagjosite 
waveform.  One  simple  application  which  baa  been  used  extenaiveiy  in  the 
past  is  that  of  determining  the  time  interval  between  the  arrivals  of 
the  Type  I  air blast -induced  wave  and  the  Type  II  oscillatory  coaponent. 

The  interval  between  the  arrivals  of  the  first  refraction  and  the  first 
reflection  can  also  be  estimated  with  fair  accuracy.  These  data  together 
with,  say,  a  triangular  representation  of  the  moving  airblast  loading 
permit  the  gross  calculation  of  the  relative  strengths  of  these  ear'y 
mot i ohs. 

If,  as  in  many  cases  it  is,  the  response  spectrum  of  the  total  motion 
of  the  ground  m  determined  by  current  ground shock  prediction  techniques 
is  provided  as  a  criterion  of  the  street h  of  the  motion  in  a  given 
application ,  the  problem  of  synthesising  acceptable  time-histories  of 
the  ground  motion  is  reduced  appreciably.  As  the  arrival  time  and  x  . 
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strength  and  fore  of  the  airblaet-in&uced  motion  can  be  estimated  with 
good  accuracy,  and  as  the  arrival  times  of  reflect ions,  refractions  and 
other  waves  as  desired  can  be  obtained  from  the  time-distance  curve,  only 
the  relative  strengths  of  these  latter  waves  must  be  estimated,  as  the 
response  spectrum  of  the  composite  wave  roust  match  the  given  spectrum. 
Obviously  a  great  deal  of  intuition  is  involved  in  assigning  even  relative 
strengths,  but  at  least  this  method  provides  a  broad  guideline  for  con¬ 
structing  the  Type  II  oscillation  which  heretofore  has  not  been  available. 
And  furthermore,  the  resulting  waveform  will  reflect  the  influence  of  the 
particular  site  and  weapon  parameters. 

The  directions  and  phasing  of  the  strong  early  motions  is  of  vital 
importance  in  the  design  of  many  hardened  structures,  in  particular  those 
extending  to  or  above  the  surface  of  the  ground.  Here,  owing  to  the  lack 
of  appreciable  stiffness  in  the  upward  direction,  the  phasing  between 
upward  and  downward  outrunning  ground  motions  and  between  these  motions 
and  the  arrival  of  the  airblast  wave  can  have  a  significant  influence  on 
the  response  of  the  structure.  Again,  while  the  strength  of  the  motions 
can  only  be  approximated  grossly,  the  range  of  phase  relationships  of  the 
first  arrivals  can  be  obtained  directly  from  the  time-distance  curve. 

It  is  clear  that  a  wavefront  diagram,  time-distance  curve  or  other  highly 
simplified  graphical  means  cannot  be  expected  to  yield  accurate,  detailed 
time-histories  of  a  phenomenon  as  complex  as  the  propagation  of  nuclear 
blast  generated  waves  in  an  inelastic,  stratified  half -space.  In  the 
absence  of  further  atmospheric  nuclear  testing,  more  accurate  solutions 
mint  corns  from  comprehensi ve  analytical  treatments  supported  by  experi¬ 
ment.  However,  until  the  results  of  such  solutions  become  available  in  a 
form  which  will  permit  the  engineer  readily  to  construct  waveforms  directly 
applicable  to  the  site  and  weapon  conditions  of  interest  to  him,  these 
graphical  techniques ,  qualitative  as  they  may  be,  offer  vitally  needed 
assistance. 
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3.  TM3CR2TICA1  AHALJfSIS  OF 


Ft  smkftfr ;  Aa:  * 


A.  Background.  Bely  waves  generated  in  an  Infinite  or  semi-infinite 
homogeneous >  elastic,  undamped  medium  by  the  action  of  a  disturbance 
applied  suddenly  at  a  point  do  not  exhibit  oscillatory  motions  unless  the 
disturbance  itself  is  oscillatory.  In  any  case,  the  input  wave  is  propa¬ 
gated  throughout  the  medium  unchanged  in  form.  In  seeking  an  explanation 
for  the  oscillations  reported  in  Reference  1,  possible  deviations  from 
these  ideal  conditions  therefore  were  investigated. 


In  Section  2,  the  possibility  was  examined  that  the  oscillatory  nature  of 
the  observed  waves  (Figure  2)  was  due  to  the  superposition  of  waves 
arriving  at  the  point  of  interest  from  a  common  source,  but  by  different 
routes.  Although  except  for  the  very  early  motions,  quantitative  corre¬ 
lations  could  not  be  obtained,  it  was  shown  that  nonuniformities  in  the  site 
properties  are  undoubtedly  responsible  to  a  large  if  not  major  extent  for 
the  oscillatory  nature  of  the  Type  II  wave.  Uonhomogeneity  of  the  site, 
however,  is  not  the  only  possible  cause  of  irregularity  in  the  waveform. 

The  energy  transmitted  to  the  ground  by  a  surface  burst  of  a  nuclear  weapon 
is  not  steady  either  in  apace  or  in  time.  The  alrblaat  wave  moving  outward 
from  the  point  of  burst  varies  in  strength  and  in  waveform,  and  as  it  is 
continuously  isparting  energy  to  a  nonlinear  medium,  deviations  freon  the 
initial  pulse  shape  will  be  generated. 


Damping  characteristics  of  the  soil  might  also  lead  to  oscillations,  for 
example,  the  theory  is  advanced  in  Reference  17  that  a  sharp  seismic  dis¬ 
turbance  in  a  homogeneous  half  * space  gives  rise  bo  a  traveling  wavelet  of  a 
shape  determined  by  the  nature  of  the  earth's  absorption  spectre*  for  elastic 
waves.  It  is  shewn  that  seismograms  can  be  duplicated  by  successions  of 
these  wavelets  generally  overlapping  but  sometimes  in  the  clear.  It  was 
further  shown  that  the  center  of  the  wavelet  travels  with  a  velocity 
characteristic  of  the  medium  and  that  the  wavelet  broadens  as  It  propagates 
outward  from  the  source.  While  this  phenomenon  may  not  be  of  principal 
significance  either  here  or  in  seismograms,  it  does  serve  to  Indicate  the 
variety  of  effects  which  may  lead  to  oscillations. 
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A  farther  possibility  it  that  the  oscillations  are  not  due  entirely  to  body 
win  but  that  they  alto  Include  action*  due  to  surface  were*  of  the  Rayleigh 
type.  Sue  atreoftha  of  the  oscillations  due  to  surface  ware*  would  be 
expected  to  be  greatest  at  the  near-surface  locations  at  which  the  type  II 
waves  were  measured.  This  possibility  was  investigated  here  and  Rayleigh 
wave  motions  resulting  from  loading  patterns  typical  of  those  generated  by 
nuclear  blasts  were  examined. 

2n  Reference  18,  previous  studies  of  magnitudes  eud  farm  of  Rayleigh  waves 
generated  by  nic lear -buret-type  loading  functions  are  reviewed.  The  complete 
solution  for  stresses  and  dlaplacamnts  produced  by  a  suddenly-applied  con¬ 
centrated  load  with  a  step  distribution  in  tine  has  been  treated  by  Pekeris 
and  IdLfaon  (19)*  The  portions  of  these  quantities  which  ar  t  due  to  Rayleigh 
effects  were  extracted  from  these  results  by  Chao,  Blelch,  and  Sackrna  (20) 
as  the  contribution  of  certain  poles  in  the  Inversion  integrals  of  the 
Feharls  solution.  Baron  and  Lectrb  (18)  extended  this  work,  considering  the 
Rayleigh  effects  due  to  a  decaying  pressure  pulse  of  uniform  strength  acting 
over  a  circular  surface  area  of  increasing  radius. 

In  this  study,  the  loading  function  wee  further  modified  to  Include  a  very 
strong  initial  pulse  Which  is  more  nearly  representative  of  the  high  yield 
nuclear  blest.  A  load  of  varying  magnitude  end  waveform  was  considered  to 
save  across  the  interface  of  an  elastic ,  homogeneous,  isotropic  half -space. 
The  ground  motion  due  to  Rayleigh  waves  alone  was  calculated  for  shallow 
depths  at  wsaricoe  ranges,  and  the  equations  for  the  body  waves  at  depths 
directly  below  the  center  of  burst  were  derived .  The  oscillatory  nature  of 
the  Rayleigh  waves  was  axial  red  sal  their  strengths  cosparad  with  the  form 
and  strength  of  the  airblart  induced  ground  shock. 

B.  Mathematical  Model. 

(l)  Aeeuaptioas.  It  was  assuaed  that  the  half-space  is  heno- 
Kumm  and  elastic.  The  pressure  dm  to  a  surface  burst  of  a  nuclear 
weapon  was  represented  as  an  instantaneous  point  loading  followed  by  a 
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radially  expanding  ring  or  circular  line  loading  of  high  intensity  with  the 
loading  within  the  ring  being  of  relatively  low  Intensity.  The  radial 
velocity  of  the  loading  front  was  assumed  to  be  Inversely  proportional  to 
the  square  root  of  the  time.  Only  the  effect  of  normal  pressure  on  the 
surface  of  the  half -space  was  considered. 

(2)  Simplified  Analysis.  The  expanding  ring  loading  or  blast 
wavefront  of  high  intensity  was  based  on  data  given  in  Reference  21  for  a 
20  megaton  weapon  and  reproduced  here  as  Figure  37  •  Here  the  radial  dis¬ 
tances  traveled  by  the  wave  and  the  peak  pressures  are  given  for  six 
locations  of  the  wavefront.  In  converting  these  data  to  a  more  amenable  form 
for  use  in  the  analysis,  they  were  first  approximated  by  a  ring  of  uniform, 
but  time-varying  distribution,  of  total  force  SQ,  and  a  ring  with  a  parabolic 
distribution,  representing  the  steep  wavefront,  of  total  force  S, 

(figure  38).  The  parabolic  ring  was  further  simplified  by  expressing  it  as 
a  Delta  function  traveling  with  a  velocity  inversely  proportional  to  the 
square  root  of  the  time.  The  expanding  circular  area  loading  was  approxi¬ 
mated  by  a  Heaviside  step  function  (Figure  39)* 

Ihe  circular  line  ring  loading  front  was  assumed  to  travel  slower  than  the 
area  loading  front  since  the  resultant  farce  of  the  formsr  is  st  s  distance 
dr/4  behind  the  leading  edge  of  the  latter. 


In  Able  13  the  data  of  Figure  37  have  been  tabulated  and  several  signifi¬ 
cant  parameters  calculated.  It  is  seen  that  the  ratio  of  S^/So  Is  fairly 
constant  for  the  range  of  primary  interest  and  waa  thus  fixed  at  0.3* 
Moreover,  the  product  PQ  -  SQd  Is  also  seen  to  be  reasonably  constant, 
averaging  about  12  x  10^5  pound  feet. 


akdSi MM 
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where  a^  is  the  wavefront  velocity  coefficient  for  the  ring  loading  wave. 
The  coefficients  a  and  a^,  based  on  the  data  of  Figure  37  are  presented  in 
Table  14. 


Table  14 

EVALUATION  OF  WAVEFRONT  VELOCITY  COEFFICIENTS  a  AND  a1  (FROM  FIGURE  37) 


t 

(ns) 

d 

(ft) 

d 

r 

(ft) 

di 

(«) 

a 

(ft2/sec  x  10^) 

*1 

(ft2/sec  x  106) 

42 

2360 

660 

2190 

33 

28 

122 

3560 

900 

3330 

26 

22 

291 

5040 

1120 

4760 

22 

19 

The  average  of  a  is  about  23  x  10^  square  feet  per  second  (Figure  40) 
while  the  value  of  a^  night  be  taken  about  20  x  ic6  square  feet  per 
second.  However,  in  view  of  the  indicated  trend  of  the  center  of  the 
ring  loading  to  lag  more  and  more  as  tine  advances ,  the  value  for  a^  used 
in  the  computer  solution  was  taken  as  only  60  percent  of  that  for  a. 

Thus, 

/ 

t 

I  a^  *  0.6  a 

As  a  first  approxlaatlon,  only  the  vertical  component  of  the  Rayleigh  wave 
was  considered  for  stations  away  from  the  loading  tone  and  near  the  surface. 
Stations  at  those  locations  are  referred  to  as  Region  I,  (Figure  4l).  For 
stations  located  directly  under  the  loading  cone  and  deep  underground,  the 
radial  horizontal  projection  r  was  assumed  to  be  snail  In  ccaparison  with 
the  depth,  z.  For  these  stations.  Region  II  (Figure  42),  the  vertical  par¬ 
ticle  velocity  is  the  quantity  of  interest. 
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figure  feO*  Comparison  of  Have  front  Velocity  of 
Figure  57  with  Simplified  value  Where 
a  «  25  x  1CF  ft2/sec. 


Figure  fcl.  location  of  Sites  ir  Region  I 
Relative  to  Load  Source. 


8* 
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Figure  1*2.  Location  of  Sites  in  teflon  H 
Relative  to  Load  Source. 


(3)  Log  of  Cosputed  Curves.  For  teflon  I,  tbs  vertical  velocity 
component  of  the  Rsylolfh  w*ve  in  jiven  if  Bpastion  3k,  (Ajpnodix  l).  *be 
dlneaslonless  velocity  coordinate  T|  is  defined: 


2*r2y 


10*  -  68.8 


2  e 


In  presentiaf  the  results  of  the  calculations,  1)  has 
the  dlaeasionlass  tine 

t  ■  ct/r 


for  the  diasnslonloss  depth  ratios 


t/T  *  {o! 


001 

010 


been  plotted  sfalnst 
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far  the  >11  tii  ionises  "spike  ratio" 


(this  ratio  to  fairly  e^astsot  throughout  the  range  covered  toy  Figure  37 
mA  liatoi  la  fSbls  13),  and  for  the  dlaaasianless  "softness"  level*.  a, 
share  a  has  be*w  interpreted  as  a  "softness"  and/or  Location  Level.  An 
iaeraaee  la  XX  Indicates  that  the  aaAlua  bacons  softer,  or  if  the  softness 
Is  unchanged  os  the  range  of  the  site  fraa  ground  zero  bee ones  Less. 

lh  the  following  tabulation,  the  figures  in  which  the  results  are  plotted 
for  serious  seines  of  sad  l  are  shown. 


*  a7^/cr 


fbr  lagloa  XX,  the  vertical  particle  velocity  coagouent  is  expressed  by 
agnation  46,  but  a—rlcal  calculations  base  not  bean  aede. 


Bard.?  Figure  45 
Bazd-Soft,  Figure  46 
Soft-Bard,  Figure  47 
Soft,  Figure  46 


I 

* 


0.001 


Bttd,  Figure  49 
Bard -Soft,  Figure  30 
Soft-Bard,  Figure  51 
Start,  Figure  52 
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C .  Discuss ion  of  Equation.  The  pameter  a  has  been  described  as  an 

Index  of  "softness"  of  the  eediua  through  which  the  wave  travels.  However, 
it  can  also  be  shown  to  indicate  whether  the  point  of  interest  lies  in  a 
•vhseisaic  or  a  superset  sale  region,  fy  definition 


sd 

a  »  er 


02) 


But 


*  -g 


tbdr 


a 


The  shear  wave  velocity  c  can  be  related  to  the  Si latat tonal  ease  velocity 

s  *  “■  “pr,“10" 


c  ■  c 


-JhL. 

P^2(l  -  v) 


end  when  Poisson's  ratio  v  is  0.25,  c  -  0.577  c  .  7  baa  been  given  pare* 

P 

vlausly  as 


>’“s(JW5/  *  l-« 


Par  a  honogensous  nsdisn,  the  transeiaalc  point  is  defines  as  tha  location 
at  which  the  blast  wave  velocity  is  equal  to  the  dllatatlonal  ease  velocity, 
that  is,  where  ▼  »  Than  at  this  point 
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fv  -  cp  \r 


In  the  cases  considered  here  responses  Sieve  been  evaluated  for  a  =  0*1; 


0  3>  0.5  »nd  1.0.  In  the  first  three  cases >  then,  the  point  of  interest 
lies  veil  beyond  the  transeisadc  point  while  in  the  last  case  the  point  of 
interest  coincides  with  the  tranteismic  point. 


Similarly ,  a  can  also  be  expressed  in  terms  of  the  Rayleigh  wave  velocity j 


cr.  Then  since  c^  *  c/y. 


and  when  v  »  c 


a)v  *  c  "  2  (I 


where  in  this  case  d  is  the  distance  to  the  point  at  which  the  velocity  of 
tho  Rayleigh  wave  coincides  with  that  of  the  eirblast  wave.  Then  the 
corresponding  time,  t^,  it 


d  ra 

5a  c 


A  staple  analysis  of  the  vertical  velocity  response  of  the  medium  at  the 
surface  will,  reveal  readily  seme  interesting  relationships  without  the  aid 
of  the  electronic  computer.  Placing  t  ■  0  in  Equation  the  following 
expression  for  the  vertical  velocity  is  obtained: 
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The  above  expression  shows  that  at  three  values  of  the  nondimensional  time 

r>  the  velocity  w  at  the  surface  approaches  infinity;  namely,  at  tq  *  y,  at 

t  *  7  -  a  ;  and  at  t  =  y  -  a,.  The  significance  of  the  three  values  of 

t  which  yield  infinite  velocities  at  the  surface  is  easily  seen  from  the 

time -distance  curve  of  Figure  43(a).  Here  two  curves  are  used  to  indicate 

the  oravel  of  the  blast  wave,  one  for  the  ring  loading  and  one  for  the 

cylindrical  loading.  The  point  on  each  curve  at  which  the  Rayleigh  wave 

velocity  Is  equal  to  the  blast  wave  velocity  was  shown  previously  to  be  at 

tijses  r  -  a  ,  and  t  ~  Q.  and  at  ranges  d/x*  =  2a  ,/y  and  d/r  *  a  h 
si  so  si  '  so' 

respectively. 


Figure  43(a) .  Tims -Distance  Curve  for  Rayleigh  Waves  Generated 
by  Cylindrical  and  Ring  Alrblast  Loadings, 
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»eoii*ansional  velocity  can 
expression: 

dMi.i 

dt  c  fit 

Then  when 


an  .  c,  iMl 

dt  #  dT 


1.0 


ikl  -  c  ,  1M1  *  i 
dt  r*  dt  7 


be  related  to  the  dimensional  velocity  by  the 


Figure  43(b),  Time-Distance  Curve 


From  the  geometry  of  figure  43(b),  it  is  seen  that 


7  Cl  "  $)p3 


where  the  subscript,  p,  indicates  values  at  the  point  of  tangency ,  p,  and 
the  subscript,  r,  indicates  values  at  d/r  =  1.0.  Then 


T  «  7  -  Q 
r 

Thus  the  Rayleigh  wave  generated  by  the  cylindrical  loading  arrives  at 
d/r  «  1.0  when 

<f  a  7  *  0 

so  '  so 

itdd  the  Rayleigh  wave  generated  by  the  ring  loading  arrives  when 

v  *  7  -  °,i 

3Eb...  relationships  are  wild  of  ':ours  only  for  ot  a  f/2.  At  lnrger  values 
of  a,  d/r  is  greater  than  1.0  and  the  point  of  tangency,  that  is  the  point 
at  which  the  airhlnst  and  Rayleigh  wave  velocitiep  are  equal,  occurs  at  a 
ground  range  beyond  the  point  of  interest. 


90 


MVL-TR-65*"S7 


It  is  seen  that  these  are  precisely  the  values  of  t  at  which  the  velocity 
w  becomes  infinite.  Thus  at  range  d  =  ?,  the  value  t  =  7  simply  narks  the 
arrival  of  a  Rayleigh  wave  which  was  generated  at  the  source  at  r  *  0. 

The  value  t  ~  7  -  C*ai  indicates  the  arrival  at  r  of  a  Rayleigh  wave  either 
which  traveled  at  the  velocity  of  the  ring  loading  until  it  became  less 
than  the  Rayleigh  wave  velocity  or  which  originated  at  the  point  where  the 
two  velocities  were  equa.I.  At  greater  ranges  where  >  v,  the  Rayleigh 
wave  continued  at  constant  velocity  c^,  arriving  at  range  r  when 
t  =  7  -  Of  . .  Similarly  the  third  infinite  value  of  v  oc cured  when  the 
Rayleigh  wave  generated  by  the  cylindrical  loading  at  the  point  where  the 
two  velocities  were  equal  reached  range  r  at  t  ■  7  -  a  . 

SO 


Figure  44.  Vertical  Velocity  at  Surface. 

Using  Equation  (13),  the  response  can  be  approximated  as  shown  in  Figure  44. 
If  the  two  pressure  loadings  were  to  move  with  the  same  speed,  the  center 
of  pressure  of  the  line  loading  would  coincide  with  the  wavefront  of  the 
area  loading  and  their  phases  would  coincide.  It  Is  improbable  that  an 


airblast  wave  having  the  continuous  variation  of  pressure  with  time  typical 

of  those  generated  by  nuclear  bursts  would  produce  the  two  distinct  pulses 

such  as  the  first  two  shown  in  Figure  44  even  though  the  pressure  decay 

rate  following  the  wavefront  were  very  rapid.  Thus  the  two  simple  peaks  at 

t  and  t  ,  would  merge  into  one. 
so  si 

Another  characteristic  of  the  re  .ponse  that  can  be  obtained  by  Equation  13 
Is  the  relatively  slow  rise  in  the  response,  w,  followed  by  a  very  abrupt 
decay.  This  phenomenon  exists  not  only  on  the  surface  hut  also  below  the 
surface  to  considerable  depth,  even  though  the  response  will  be  finite  due 
to  the  presence  of  the  tern  t  (7^  -  l/5)^^  in  Equation  34.  The  presence 
of  this  spatial  decay  may  be  of  importance  in  regard  to  spalling  of  rock. 

D.  Calculated  Responses,  Region  I.  Responses  obtained  by  computer 
solutions  of  Equation  34  are  shown  in  Figures  45  through  52.  It  may  be 
noted  that  in  all  eight  erses,  the  form  of  the  vertical  velocities  due  to 
Rayleigh  waves  agrees  closely  with  that  predicted  in  Figure  44  and  by  a 

study  of  the  relationship  of  the  parameters.  In  each  case  a  wave  origina- 

\ 

ting  at  the  source  at  *f  »  0  reaches  range  d  •  r  at  t  ■  7  0  1.08  and  at 
t  ■  J  -  a  ,  o  negative  spike  of  lower  amplitude  occurs.  This  latter  wave 
originates  at  the  rang*.;  at  which  the  velocity  of  the  airblast  wave  from  the 

ring  loading  equals  the  Rayleigh  wave  velocity,  i.e.,  at  d/r  *  20^/?. 

\ 

The  first  wave  to  reach  the  point  of  interest  in  all  figures  except  48  and 
52  reached  range  d-r  at  r  -  7  ■  «g0  *nd  originated  at  range  d/r  -  20^/7 

at  <r  m  Ct  when  the  velocity  of  the  airblast  wave  from  the  cylindrical 

so 

loading  was  equal  to  the  Rayleigh  wave  velocity.  In  Figures  48  and  52, 

the  values  of  3  and  a  .  ere  1.0  and  0.6  respectively.  Thus  Rayleigh 
80  Si 

92  \ 


-  ~i  n>l—  f  l«  k  m 


so 


AFWL-TR-65-67 


"... 

'•"'1 
V  **j»  * 


c  OT'O 


Rayleigh  Waves  in  Figure  50.  7)  vs.  t  for  Rayleigh  Waves  in 
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v»ret  arrived  at  d  »  r  simultaneously  with  the  airblast  wave*.  Tbeae  tlm»g 
were, 

To  *  Jo”  *  °‘29 
80 

for  the  cylindrical  loading  and 

.2 

T1  =  BT  '  oM 

sl 

for  the  ring  loading. 

It  nay  be  noted  that  the  Rayleigh  waves  resulting  frees  the  passage  of  the 
cylindrical  loading  across  range  r  are  not  shown  in  Figures  ^  anl  52  as 
they  arrive  prior  to  the  niniaun  value  of  T  shown. 

In  all  these  curves,  based  on  an  elastic  madiua,  a  fixed  spike  ratio 
Si/S0  °  0.3  and  e  fixed  velocity  parameter  ratio  *  0.6,  there 

should  be  no  fundamental  changes  in  shape  if  the  ratios  are  varied  only 
slfghtly.  As  shown  previously,  the  phasic*  between  the  Rayleigh  waves 
resulting  from  the  cylindrical  and  from  the  rii*  lowdii*  is 

A*  "  “sl  ‘  “so 

while  the  phasing  between  the  Rayleigh  wave  due  to  the  ring  loading  and 
that  ge .aerated  at  the  source  is 

4’  '  V 

Both  of  theoe  difference,  ere  vnlld,  of  couree,  only  if  both  a  nod  a 
,  so  sl 

are  equal  to  or  greater  than  y/2, 

Tims -distant:;  carves  for  each  of  the  valuta  of  are  shown  in  Figures 
53  through  56.  Also  indicated  in  the  figures  are  the  arrivals  of  the 
di  la  tail  octal  waves. 
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Figure  55.  Time-Distance  Curve  for  Rayleigh  and 

Dilatational  Waves,  a  *  0.50. 
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Figure  56.  Time-Distance  Curve  for  Rayleigh  c  d 

Dilatations!  Waves ,  a  »  1,00. 
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Vmm  *jeak#  would  be  attenuated  significantly  due  to  two  factors.  First# 
in  the  region  surrounding  the  point  of  burst,  the  elastic  theory  is  invalid 

since  at  the  extremely  high  pressures  in  this  region  the  medium  acts 
essentially  like  a  fluid,  in  which  of  course,  Rayleigh  waves  do  not  occur. 
Thus  waves  originating  within  this  region  would  be  expected  to  be  less 
severe  than  those  predicted  by  elastic  theory. 


A  second  factor  tending  to  reduce  the  strength  of  the  waves  is  that  even 
cce§>®tant  rock  will  not  remain  elastic  at  the  high  stresses  accompanying 
th<  high  velocities.  Soils,  of  course,  will  produce  even  greater  atten¬ 
uations  particularly  if  near  the  surface,  since  they  will  yield  appreciably 
in  compression  and  can  sustain  tensile  loads  no  greater  than  the  in  situ 
ccBgtr&*slon. 


Thus  the  strengths  of  the  waves  shown  in  the  curves  represent  an  upper 
bound  of  those  that  might  be  expected  to  occur  in  a  real  medium.  Further, 
the  wave  velocities  and  shapes  would  be  altered  by  nonlinearities  in  the 
medium,  thus  changing  the  arrival  times  and  phasing  of  the  waves  from  those 
indicated.  This  latter  effect  however  would  be  expected  to  be  relatively 


small. 


Although  the  strengths  of  Rayleigh  waves  in  a  real  medium  could  be  appre¬ 
ciably  less  than  those  shown  in  the  curves,  it  is  of  interest  to  compare 
their  strengths  with  those  of  dilatational  waves  generated  by  the  airblaet 
loading  at  the  same  sites.  Of  the  four  cases  far  which  numerical  values 


were  obtained,  one  (a^  »  1.00)  was  at  the  transeismie  point  while  the 


remaining  three  were  in  the  subseismlc  or  outrunning  region.  In  the  tran- 
selsmic  and  subseismlc  regions,  the  vertical  velocity  at  the  surface  due  to 
an  alrblaat  pressure  p  is  not  equal  to  p/pc  but  is  Increased  by  the  con¬ 
tribution  of  the  ground-transmitted  shock  from  the  eirblast  wave  when  It 
was  at  a  lesser  range.  Reference  26  recommends  the  use  of  a  factor  of  1.4 


to  account  for  this  effect.  Thus 


where  u  Is  the  vertical  velocity  due  to  the  airblast  pressure 
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In  calculating  the  strengths  of  the  Rayleigh  waves*  it  was  assumed 


where  8  was  the  total  force  exerted  on  the  surface  of  the  ground  fey  the 

V 

cylindrical  loading  and  p  was  the  overpressure  due  to  the  cylindrical 

w 

loading.  Referring  to  liable  13*  it  is  seen  that  the  pressure  due  to  the 
ring  loading  is  about-  three  times  p  ,  Then 

O 


Letting  p  -  p.,  and  d  *  r,  c  ~  0-577  c 
—  P 

k.2  SQ  2.k  SQ 

i  -  - -  =  — 

«r  pCp  «r  pc 


The  ground  velocity  parameter  1)  due  to  the  Rayleigh  wave  was  given  as 

7)  „ 

cS 

o 

Substituting  ep  »  u/c 

•n  62.8ocrS> 

1  So 


Then 

*  JL 

i  ■  it 

Now  in  Figures  L  5  through  52  it  is  seen  that  the  peak  strengths  of  the 
liayleigh  waves  generated  by  the  ring  loading  are  very  high*  exceeding 
200  units.  Thus  for  the  case  of  an  elastic  medium*  and  recognizing  cha 
approximate  nature  of  the  comparison*  It  lb  seen  that  the  strength  of  the 
Rayleigh  wave  due  to  the  ring  loading  can  be  larger  than  that  due  to  the 
dilatations!  wave  generated  by  the  ring  loading. 
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E.  Rayleigh  Waves  la  Real  Media.  A  direct  quantitative  comparison 
of  the  Say Leigh  wave  pattern  resulting  from  a  disturbance  in  an  elastic, 
homogeneous  medium  with  the  oscillatory  motions  observed  at  highly  strati¬ 
fied  real  sites,  of  course,  is  not  possible.  This  work  however  was  directed 
primarily  toward  the  identification  of  phenomena  which  might  contribute 
significantly  to  the  Type  II  oscillation  rather  than  toward  direct  corre¬ 
lations.  In  this  respect  it  appears  that  the  investigation  was  successful 
and  that  several  points  of  importance  can  be  established. 

The  Rayleigh  wave  generated  by  a  step  loading  consists  of  a  single  pulse¬ 
like  motion  in  the  vertical  direction  with  negligible  oscillation.  The 
wave  motion  resulting  from  an  impulse  loading  appears  to  consist  of  a 
single-cycle  oscillation.  In  neither  case  is  a  train  of  waves  generated. 

It  is  doubtful  if  the  wave  transmitted  directly  from  the  initial  distur¬ 
bance  would  be  of  significance  in  a  real  medium.  As  it  must  follow  the 
direct-induced  congressional  wave  and  the  air-induced  wave,  it  must  pass 
through  t;.e  ruptured  zone  surrounding  the  crater.  In  this  region  it  is 
likely  that  the  highly  inelastic  behavior  of  the  medium  would  damp  the  wave 
to  a  negligible  strength. 

Even  for  the  elastic  case  it  may  be  noted  that  the  wave  generated  by  the 
steep-fronted  airblast  wave  as  it  recedes  from  ground  zero  at  a  velocity 
equal  to  that  of  the  Reyleigh  wave  is  by  far  the  most  severe.  As  the 
velocities  of  the  airblast  wave  and  the  Rayleigh  wave  are  matched  outside 
of  the  rupture  zone,  and  as  this  Rayleigh  wave  is  generated  in  a  medium 
which  has  been  preconsolidated  by  the  passage  of  the  iilatational  wave,  it 
is  possible  that  the  strength  of  this  wave  as  predicted  by  elastic  theory 
may  not  be  too  far  in  error. 

The  form  of  an  actual  airblast  wave  consists  of  a  single  exponentially 
decaying  pulse  rather  than  the  two  distinct  blast  loadings  assumed  in  the 
analysis.  Thus  the  two  Rayleigh  waves  generated  by  the  blast  loadings 
would,  in  fact,  be  superinposed .  In  addition,  the  strength  of  the  combined 
wave  would  be  less  than  that  shown  due  to  dissipation,  the  oscillation  would 
be  brc&lened  due  to  the  nonlinearity  of  the  medium  and  the  symmetry  of  the 
oscillation  would  be  distorted  by  the  inability  of  the  medium  to  sustain 
large  tensile  stresses. 
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In  this  analysis ,  it  was  assumed  that  the  medium  was  homogeneous.  At 
layered  sites,  however,  and  for  certain  ratios  of  density  and  Lame’s  cons¬ 
tants  of  the  two  media,  a  modified  form  of  the  Rayleigh  wave  can  exist  at 
the  interface.  As  it  travels  at  a  velocity  slightly  different  from  that 
of  the  surface  Rayleigh  wave,  it  would  be  expected,  that  the  motion  at 
shallow  depths  would  be  the  superposition  of  the  motions  due  to  the  two 
sets  of  Rayleigh  waves  but  at  some  difference  in  phase.  Layering  would 
also  be  expected  to  result  in  other  phenomena,  such  as  refractions  and 
reflections  at  the  interface,  which  were  not  considered  here. 

Quantitatively,  then,  a  rough  picture  of  the  Rayleigh  waves  at  a  real  site 
can  be  constructed.  In  the  outrunning  or  subseismic  region  where  the 
Type  II  oscillations  were  observed,  the  Rayleigh  wave  can  never  arrive 
first,  although  near  the  transeismic  point  in  a  homogeneous  medium  it  may 
arrive  simultaneously  with  the  airblast-induced  motion.  The  initial 
direction  of  the  vertical  conponent  will  be  upward  similar  to  the  initial 
upward  motion  of  the  Type  II  wave.  The  principal  part  of  the  ground 
velocity  due  to  the  Rayleigh  wave  will  consist  of  one  full  cycle  oscil¬ 
lation. 

The  lack  of  an  indication  of  the  intensity  of  the  Rayleigh  wave  in  real, 
layered  media  precludes  the  positive  identification  of  any  contribution 
they  may  have  made  to  the  observed  Type  II  waveforms.  Since  they  arrived 
at  the  point  of  measurement  after  the  arrival  of  the  refracted  wave,  they 
could  not  have  been  responsible  for  the  initial  upward  motion.  It  is 
possible  that  some  of  the  later  oscillations  however  are  due  at  least  in 
part  to  Rayleigh  waves. 
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APPENDIX  A 

CONSTRUCTION  OF  A  REFRACTION  WAVEFRONT  DIAGRAM 

The  refraction  wavefront  diagram  is  a  graphical  representation  in  a 
vertical  plane  passing  through  the  source  and  the  station  of  interest  of 
the  positions  at  discrete  tine  intervals  of  the  direct-induced  and  refrac¬ 
ted  waves  emanating  from  the  source.  Since  these  waves  will  always  precede 
reflected  or  shear  waves,  the  refraction  wavefront  diagram  indicates 
directly  cho  earliest  motion  at  any  point.  With  the  assumption  of  linearity 
of  the  medium,  later  arrivals  can  be  depicted  in  a  similar  manner  as  can  the 
propagation  of  reflected  waves. 

A  great  degree  of  detail  can  be  included  in  the  wavefront  diagram  if  it  is 
desired  to  do  so.  Nonhomcgeneity  can  be  represented  by  many  homogeneous 
layers  of  very  small  thickness,  the  layers  may  be  oriented  in  any  manner 
desired,  and  inclusions  of  any  shape  or  "site  may  be  considered .  r or  the 
refraction  diagram,  the  assumption  of  linearity  of  the  medium  is  not 
necessary  as  the  first  waves  to  arrive  at  a  point  will  always  travel  at  the 
seismic  velocity  of  the  medium.  For  calculating  later  arrivals  through 
nonlinear  media,  certain  rules  of  thumb  have  been  employed,  although  it  is 
in  this  respect  that  the  most  severe  limitation  of  the  wavefront  diagram  is 
evident;  the  strength  of  the  wave  cannot  be  detenained. 

The  practical  limit  to  the  detail  that  should  be  included  in  the  wavefront 
diagram  is  dependent,  of  course,  on  the  accuracies  with  which  the  site 
profile,  the  soil  properties,  the  weapon  size  sod  coupling  efficiency  and 
the  burst  location  are  known.  Since  in  most  cases  these  parameters  are 
known  only  within  broad  ranges,  many  refinements  are  not  warranted. 

In  the  example  d!«-us«ed  here,  a  three-layered  system  is  considered 
(Figure  58)-  The  upper  or  surface  layer,  with  seismic  velocity  ■  5000 
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feet  per  second,  extend*  to  *  depth  of  2500  feet;  the  second  'v«)  ’-ith 
Vg  »  7500  feet  per  second  is  UOOQ  feet  in  thickness;  suid  the  third  o 
bottom  layer  ha*  a  seismic  velocity  Vj  ■  10,000  feet  per  second, 
the  two  Interfaces  are  assumed  to  lie  parallel  to  the  surface  r.V.^agh,  as 
noted  above ,  this  as  suction  is  not  necessary  and  any  o.rie-;  .r  ;...a  of  the 
layers  could  be  considered. 

STEP  1:  The  scale  is  first  selected  and  the  soil  profile  drawn  in.  To 
facilitate  construction,  the  horizontal  scale  (range)  should  be  equal  to 
the  vertical  scale  (depth).  Yet  at  most  sites,  the  range  of  interest  ray 
extend  to  several  thousands  of  feet  while  detail  may  be  desired  in  a 
relatively  thin  layer.  Moreover,  the  diagram  construction  begins  at  the 
source  and  precedes  outward  thus  accumulating  erroru.  A  fairly  large 
scale  then  is  essential  if  any  accuracy  is  to  be  achieved. 

Beyond  the  range  at  which  a  refraction  from  the  deepest  layer  returns  to 
the  surface  the  shape  of  the  wavefront  and  its  velocity  remain  constant. 
Thus  the  wavefront  diagram  need  be  extended  only  to  this  range  which  nay  be 
caeqputed  by  Equation  8  or  derived  from  the  relationships  given  in  Section  2. 
That  la 


and 


for  a  two-layered  system 


x 

c 


for  a  three-layered  system 


The  relationship  for  the  two-layered  system  is  shown  in  Figure  57. 
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ligurc  57.  Range  -  Thickness  Ratio  at  which  a  Refracted  Wave 
First  Reaches  the  Surface  in  a  Two-Layer  System. 


STEP  2:  Rays  at  the  critical  angle  beyond  which  all  waves  are  reflected 
at  the  interfaces  are  then  drawn.  In  Figure  58  the  rsy  OA  represents  the 
limiting  path  in  the  Layer  1  -  Layer  2  system.  The  critical  angles  are 
determined  according  to  Snell's  Law  as  described  in  Section  2. 

STEP  3:  Select  a  tine  Interval  uwd  from  the  source  "0”,  draw  in  Layer  1 
circular  arcs  representing  the  position  of  the  direct -induced  wave  at  these 
times.  The  arcs  should  be  extended  lightly  into  Layer  2  as  shown  by  the 
dashed  lines  in  Figure  58. 


STEP 


The  wavefronts  in  Layer •»  2  and  3  will  not  be  circular  arcs. 


Directly  below  the  source  the  wave  will  propagate  vertically  at  the  seismic 
velocity  of  the  layer.  Along  th  interface  with  a  shallower,  softer  layer 
and  at  angles  from  the  source  less  than  the  critical  angle,  the  horizontal 
component  of  the  wave  velocity  ir.  the  lover  Liyer  will  be  less  than  the 
horizontal  component  of  the  velocity  in  the  upper  layer.  At  Points  A  and 
B,  the  horizontal  v c  i  or  1  % y  c  tY ^ nts  ip,  the  two  adjacent  layers  will  oe 


exactly  equal  and  at  greater  ranges ,  the  velocity  in  the  lower  layer  along 
the  interface  will  exceed  that  in  the  upper  layer.  Tiius  along  an  interface , 
between  Layers  1  and  2  for  example,  and  at  ranges  greater  than  Point  A,  the 
initial  motion  is  generated  by  the  first  disturbance  reaching  Point  A. 
However,  at  all  points  within  Layer  2  below  the  interface  with  Layer  1, 
except  for  those  regions  where  waves  are  refracted  into  Layer  2  from  lower 
layers ,  the  initial  motion  is  generated  by  disturbances  arriving  at  the 
interface  from  ranges  less  than  Point  A.  As  noted  earlier,  then,  the  wave 
profiles  in  Layer  2  are  not  circular  arcs. 

At  ranges  shorter  than  those  indicated  by  the  critical  ray,  intermediate 
points  on  the  wavefront  lying  between  the  interface  and  the  depth  scale 

can  be  determined  graphically  by  constructing  radii  from  the  Bource  through 
the  region  of  interest,  such  as  line  abc.  As  the  dashed  line  is  simply  an 
extension  of  the  circular  arcs  in  Layer  1,  3b  is  the  distance  the  wave 
would  travel  in  Layer  2  if  Vg  were  equal  to  V^.  Then  the  actual  distance 
abc  must  he  ah  x  Vg/\L  .  Three  or  four  such  radii  are  sufficient  for  the 
-.c curacy  needed  in  most  practical  applications. 

At  ranges  greater  than  Point  A  for  Layer  2  and  Point  B  for  Layer  3,  the 
wavefronts  V  -.  u .  euiuf.ir4-:  .  L.  . 

A  and  B  respectively. 

STEP  3:  Hie  fronts  in  layer  1  along  the  interface  between  Layers  1  and  2 
beyond  Point  A  are  first  reached  and  set  in  motion  by  the  disturbance 
traveling  in  and  along  the  top  of  Layer  2.  This  "underlayer  wave"  has  a 
straight  wavefront  in  layer  1  with  a  elcpe  of  V^/Vg  as  shown  in  Figure  58. 

Hie  intersection  of  the  underlayer  wave  and  the  direct  wave  having  the 
corresponding  lime  is  a  point  on  the  coincident-time  curve  AC  in  Layer  1. 

Hie  coincident-time  curve  originates  at  Point  A  in  Layer  1,  and  at  Point  B 
in  Layer  2.  In  Layer  2,  the  coincident-time  curve  represents  the  inter¬ 
section  between  the  refracted  wavefront  of  .aycr  2  a-i  the  underlayer  wave 
originating  f~oa  layer  ?, 
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The  points  where  the  coincident -time  curve  intersects  the  ground  surface 
is  represented  on  the  time-distance  curve  by  an  abrupt  change  in  the  slope 
of  the  line,  as  shown  by  Points  C  and  D  in  Figuro  8. 


The  wavefront  diagram  started  in  Figure  58  if  shown  in  its  completed  form 
in  Figure  5. 


x 
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APPENDIX  B 


THEORETICAL  ANALYSIS  OF  GROUNDSHOCK  WAVEFORMS 


Blast  Wave  of  a  Nuclear  Detonation  and  its  Idealized  Approximation 

The  characteristics  of  the  blast  wave  due  to  a  nuclear  detonation  can  be 
shown  roughly  as  in  Figure  37  (21 'L  The  maximum  pressure  occurs  at  the 
cylindrical  wavefront  from  which  it  decays  rapidly  and  levels  off  toward 
the  epicenter.  If  S  denotes  the  total  load  on  the  ground  surface  exerted 
by  the  blast  wave,  d,  the  distance  of  the  wavefront  from  the  epicenter,  r, 
the  radial  distance  of  the  site  from  the  epicenter  and  if  the  peak  pressure 
spike  at  the  wavefront  is  approximated  by  a  6- function,  then  the  normal 
pressure  a  exerted  on  the  ground  by  the  blast  wave  may  be  written  approxi¬ 
mately  as: 


<T  “  -  S, 


H  (d  -  r 


(U»> 


h(x),  B(x)  represent  the  Heaviside  step  frncti  u  'd  -  function 
respectively, 


S  +  S. 
o  1 


(l1*) 


So  "  f 


The  velocity,  v,  of  the  wavefront  decays  with  time  and  eventually  approaches 
zero  as  limit.  Because  of  the  numerical  values  given  in  Figure  37#  it  is 
assumed  to  he  inversely  proportional  to  the  square  root  of  time  counted  from 
the  lnetant  of  detonation.  Namely, 
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1 


i 


0-5) 


vhere  aim  proportional  constant.  Accordingly*  the  location  of  the 
wavefront  at  time  t  can  be  easily  computed  as: 


d 


«  2  v/at" 


(16) 


It  may  be  noted  that  this  result*  d  ~t0,5>  is  very  close  to  G.  I.  Taylor  s 
intense  blast  wave  theory  (22)  which  shews  d  ~  t0,  .  A  conqparison  between 
Equation  l6  end  the  observed  values  from  Figure  37  are  plotted  in  Figure  40. 
It  my  be  seen  that  the  agreement  is  quite  good  if  "a"  is  taken  to  be 
25  (kilofeet)  /second. 


Substitution  of  Equation  l6  into  Equation  l4  .gives  the  loading  function 
due  to  a  nuclear  blast  wave  as: 


(if) 


Mathematical  Analysis 


Assume  that,  the  ground  is  a  homogeneous*  elastic  half-space  of  density  p* 
shear  modulus  ^  Poisson's  ratio  v  -  0.25.  Using  cylindrical  coordi¬ 
nates  ac  shown  in  Figure  59*  letting  q  and  w  denote  the  radial  and  the 
vertical  displacements  respectively  and  following  las to  (23)*  these  dis¬ 
placements  may  be  derived  from  a  dilatational  potential  0  and  from  u  shear 
potential  ♦  according  to 


q 


j* 

Sr 


it. 


(10a) 
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w 


4 

dr 


1 

r 


SL 

ar 


(1Gb) 


The  governing  differential  equations  of  the  medium  nay  be  written  in  terms 
of  0  and  ♦  as: 

c  Vo  -  h  (19») 

p  at4 

=Vt  -  i\  d9b) 


* 


Figure  59.  Cylindrical  Coordinates. 

where  , 

c  •  (ji/d)1/2  and  cp  •  •  JTc  (20) 

denote  the  eh-*;*  and  diUtatiooal  wave  speeds  of  the  sMdiua  respectively. 
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The  solutions  to  Equations  19a  and  1<&  are  subjected  to  the  boundary 
conditions 

o  t  >  0,  z  «  0  (21a) 

0  t  >  0,  z  -  0  (21b) 

where  <r  is  defined  in  Equation  17  and  a  ,  <x  sre  stress  conponents  of  the 

zz  zr 

medium.  The  displacement  and  therefore  the  potentials  also  must  vanish  at 
infinity.  The  initial  conditions  are  such  that  all  quantities  must  be  zero 
for  t  -  0. 


Let  q,  w,  0,  ?,  a  and  Q,  W,  #,  f,  £  denote  Laplace  and  Laplace-Hankel 
transforas  of  q,  v,  0,  f,  or  respectively  according  to 


/*  • 

(q*  Vi  o)  m  f  (q»  v*  ff»  t, 

rft 


r)-Pt 


dt 


(q*  v,  0,  ♦,  »  /  (q,  v»  ?,  ¥, 


(22a) 


(22b) 

(23a) 

(25b) 


where  p  and  f  are  Laplace  and  Bankel  transform  parameter*  and  B  indicates 
the  path  of  integration  for  inverse  Laplace  transform  such  that  q,  v,  ... 
are  analytical  functions  of  p  to  the  right  of  Br. 

The  transformed  equations  corresponding  to  Equations  19a,  19b,  21a,  and 
21b  are: 

118 


AFWL-TR-65-67 


where 


|  j  - 

0 

(2!*a) 

|  y  * 

0 

(2lb) 

(2^  +  ks)  *  * 

2^ 

g  -  2/ji*  *  -  0 

(25a) 

2  H  + 

♦  k: 

")  y  ■  0#  z  ■  0 

(25b) 

b  •  p/c  ■  • 
P 

— JL_ 
v'ST 

;  k  •  p/c 

(26) 

and  £  is  the  lap  lace -Henkel  transform  of  the  loading  function  a  which 
takes  the  fora  (2^)  of: 


E  («»  P)  « 


!k-  .-^/P 


(27) 


# 


With  the  aid  of  the  vanishing  boundary  conditions  at  infinity#  f  and  y  can 
be  obtained  explicitly.  Substituting  into  the  transformed  ver 'too  of 
Equations  18a  and  18b#  the  transformed  vertical  and  radial  displacements 
of  the  medium  are  given  by: 


•  «pj-*  ♦  b2jl/2J-  «xp|-t  (f  *  k2jl/2|  (28) 


W  ({#p,t) 
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Where 

M  (?,P)  -  kt?(f  *  h^2  ^  *  K^/2  (5C) 

end  the  branch  of  square  root  is  chosen  to  yield  a  positive  real  part  in 
order  to  insure  the  vanishing  displacement  conditions  at  infinity. 

Front  Equations  28  and  29,  similar  transformed  expressions  for  radial  and 
vertical  velocities#  stress  and  strain  components  can  be  easily  derived. 

In  particular,  the  most  lsg>ortant  quantity  needed  for  detection  of  the 
Type  II  waveform  is  the  vertical  component  v  of  the  particle  velocity  of 
the  medium,  whose  lap  lace -Hankel  transform  can  be  written  down  readily  as 
P 

Formal  Inversion  is  accomplished  by  the  aid  of  Equations  22b  and  2Jb. 

The  vertical  particle  velocity  of  the  medium  then  Is  given  by: 

*  '<r,t'i> '  MX  5P  V({,p,t)  •**  ($r)  dp  df  (Jl) 

r 

Complete  solutions  of  this  integral  valid  for  all  stations  in  the  medium 
require  extremely  lengthy  numerical  confutations.  The  present  investiga¬ 
tion  therefore  is  confined  to  two  regions  of  practical  interest  for  which 
the  ratio  s/r  is  either  very  small  or  very  large.  Because  of  the  extremal 
values  taken  for  this  ratio,  approximations  based  on  physical  argument  or 
mathematical  simplifications  are  esgloyed  to  obtain  solutions  which  can  be 
computed  numerically. 
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Solution  for  Shallow  Depth  of  Station  -  Region  I.  z/r  «  1 


For  stations  far  fron  the  loading  zone  and  near  the  surface,  see 
Figure  ll>  the  major  response  of  the  medium  is  due  to  Rayleign  waves. 

This  has  been  discovered  by  Rayleigh  and  is  well  confirmed  by  sei sinolog leal 
experience.  It  has  al30  been  found  (19)  that  if  the  transform  approach 
were  used  (as  the  present  analysis  did),  the  Rayleigh  wave  effect  can  be 
obtained  from  the  transformed  solution  by  utilizing  only  the  residues  of 
the  inversion  integrals  at  certain  poles  (Rayleigh  poles),  ignoring  branch 
integrals  which  would  occur  in  a  complete  solution. 


The  Rayleigh  poles  of  the  transformed  vertical  velocity.  Equations  20  and 
31,  are  the  zeros  of  Equation  30,  M(5,p)  -  0,  which  may  be  easily  recognized 
as  Rayleigh's  frequency  equation  of  an  elastic  half -space.  The  location  of 
the  poles  are 

p  «  ic 5/7  (32a) 

where 

7  -  |(3  *&f/2  (52b) 

ft*  contributions  from  these  two  poles  can  be  computed  with  the  aid  of 
simple  residue  theory  as 
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(34«) 


(3  4f) 


Equation  34a,  the  Rayleigh  phi»*u.  of  the  verticr.*  velocity  of  the  medium, 
is  plot  .<t  1  in  Figures  45  to  r2  i  i  w*  *ch  tfc&  ordinate  is  the  nondiaenaional 
vertical  velocity  and  the  abncijsa  is  the  aondiasnaional  time  r*  The 
meanings  of  various  p«iweters  in  +hese  J,  ratings  are  explained  in  Figures 
4l  and  42. 


Solution  for  Large  Pep H;  of  Statj on  - 


L*  z/r  »  1 


For  stations  near  the  loading  r.cne  and  deep  underground,  r ,  the  radial 
distance,  may  he  assumed  to  be  sero  in  comparison  with  the  depth,  see 
Figure  42.  The  Laplace  transform  of  the  vertical  displacement  then  nay  be 


written  as 


iir  (o,p,s)  ■  jf  "  w  (o,t,s)  e-pt4t 


/  *  W  (5,p,r)  5d5 


according  to  Equations  22a  and  23b  with  (§r)  evaluated  at  r  «  0.  Let  a 
change  of  variable 


-  22E  -  kx 

c 


he  Introduced  in  Equation  33.  Substitution  of  Equations  27  and  28  yields 


*  ««¥Mi 


*  1  1 
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w 


where 


and 


Recalling  that  w  is  the  Laplace  transform  of  the  vertical  displacement  v 
Equation  22a,  It  is  seen  that  the  quantities  1^  Ig,  in  Bqtiations  UOa, 
Uot,  therefore,  are  the  Laplace  transforms  of  some  function,  say 
Kg  such  that 


Ix  (p>*»P)> 


(P»z 


-]  i 


Kx  (t,z,P),  Kg  (t,t,p)  e 


*P)]« 


-pt 


dt 


(M) 


It  is  therffcre  apparent  that  the  detarminat i on  of  the  vertical  valtcity  d 
hinges  on  thi  determination  of  Kx  and  Kg  of  Equation  kl.  To  this  end  let 
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is  a  nundimensional  tine.  It  may  be  noted  that  the  form  of  Equation  teh 
is  expressed  in  such  a  way  that  no  ambiguity  will  arise  when  P  -  0. 
Furthermore,  all  square  roots  are  chosen  to  be  positive. 


Substitution  of  Equations  tea,  teb  and  tec  into  Equation  tea,  with  some 
mathematical  manipulation  yields  the  result. 


186 


Am-TR-65-67 


Comparing  this  expression  with  that  of  Equation  Hi,  it  can  be  seen  that 


Kx  (t,z,p)  -  0 


for 


T1  > 


1 

VT 


where  *.  and  are  defined  in  Equations  H2b  and  H2c 
Similarly,  let 


+  1  +  0x2 


t 


ik 


in  Equation  HOb  which  leads  to 


Am-TJt-6S-67 


Kg  (t,z,P) 


4T]P  +  1  +  4P2  +  ^fcTlP  +  1  +  4P2 


(46) 


for  ^  >  1 


Xg  (t,z,P)  -  0 


for  Tl  <  1 


where  x  is  defined  in  Equation  45b  * 

Referring  back  to  Equations  22a j  59>  40a,  and  4Cfb,  it  can  be  seen  that  the 
vertical  displaceaent  of  the  aedlun  directly  below  the  loading  zone  is 

*<0.t,«)  .  -  jij  |s0 

♦  Sl['Cl  it'*'°1)-,cs  (t'l'ai)]|  W) 

Hhere  K^»  Kg  are  defined  in  Equations  44  and  46  respectively.  The 
vertical  velocity  of  the  aediua  can  be  easily  obtained  by  a  direct 
differentiation*  neatly 

v  (o,t#a)  •  ^  v  (o,t*s)  (48) 

m 
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11.  ABSTRACT  ~~ . . 

A  recent  review  of  nuclear  groundshock  test  data  revealed  oscillatory  motions  of 
undetermined  origin,  but  of  sufficient  magnitude  to  be  of  concern  in  the  design  of 
shock  isolation  systems  for  underground  protective  structures.  Response  spectra 
of  the  oscillations  exhibited  amplification  ratios  far  exceeding  those  employed  in 
most  groundshock  prediction  methods,  the  ratios  being  functions  of  the  number  of 
oscillations,  and  the  amplitude  and  period  of  each  cycle.  As  the  source  of  the 
oscillations  was  not  identified,  there  was  no  rational  basis  for  relating  these 
waveform  parameters  to  such  basic  site  and  weapon  conditions  as  yield*  range, 
geologic  structure  and  properties  of  the  medium.  In  this  investigation,  two 
possible  sources  of  the  oscillatory  motion  are  examined.  First,  the  propagation 
of  waves  in  a  stratified  site  are  studied  and  their  directions  and  phase  relation¬ 
ships  estimated  by  uae  of  wavefront  diagrams  and  time-distance  curves.  Second, 
the  form  and  strength  of  Rayleigh  waves  in  an  elastic,  homogeneous  half-space 
which  result  from  surface  pressure  distributions  similar  to  those  generated  by 
nuclear  bursts  were  calculated.  In  both  cum,  oscillatory  phenomena  can  be 
predicted  and  certain  features  related  to  the  observed  oscillations.  However, 
the  simple  approaches  employed  in  this  analysis  will  not  yield  realistic  wsva 
strengths  and  thus,  the  composite  waveform  at  a  point  in  tha  half-space  cannot  be 
determined  quantitatively. 
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